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Abstract. We study the partition function for the three-colour model with domain wall 
boundary conditions. We express it in terms of certain special polynomials, which can be 
constructed recursively. Our method generalizes Kuperberg's proof of the alternating sign 
matrix theorem, replacing the six-vertex model used by Kuperberg with the eight- vertex-solid- 
on-solid model. As applications, we obtain some combinatorial results on three-colourings. 
We also conjecture an explicit formula for the free energy of the model. 



1. Introduction 

An alternating sign matrix is a square matrix with entries 1, —1 and 0, such 
that the non-zero entries in each row and column alternate in sign and add up 
to 1. Mills, Robbins and Rumsey jMRR] conjectured that the number of n x n 
alternating sign matrices is 



l!4!7!---(3n-2)! 
" n!(n + l)!---(2n-l)!' ^ ' ' 



This conjecture was proved by Zeilberger |Zlj . Soon afterwards, a much sim- 
pler proof was found by Kuperberg p<], using the six-vertex model of statistical 
mechanics. 

The six-vertex model is an example of an ice model, whose states can be iden- 
tified with what we call ice graphs; see (12. 2 p below. Alternating sign matrices 
can be identified with ice graphs satisfying domain wall boundary conditions. For 
the six-vertex model, there is a closed formula for the corresponding partition 
function, the Izergin-Korepin determinant [U lICK] . Kuperberg observed that in 
a special case, when all parameters of the model are cubic roots of unity, the par- 
tition function simply counts the number of states. He could then prove (11.11) by 
computing the corresponding limit of the determinant. 

There is a natural two-parameter extension of the six-vertex model known as 
the eight-vertex-solid-on-solid (8VS0S) model. This model was introduced by 
Baxter |B2] as a tool for solving the eight-vertex model. It is elliptic, that is. 
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the Boltzmann weights are eUiptic functions of the parameters. We stress that, 
in contrast to the eight-vertex model, the 8VS0S model is an ice model. In 
particular, for domain wall boundary conditions, its states can be identified with 
alternating sign matrices. 

It is natural to ask what happens to the 8VS0S model under Kuperberg's 
specialization of the parameters. The answer turns out to be very satisfactory: 
it degenerates to the three-colour model. It is an observation of Lenard that 
ice graphs are in bijection with three-colourings of a square lattice, such that 
adjacent squares have distinct colour [L]. The three-colour model is defined by 
assigning independent weights to the three colours; the partition function is simply 
the corresponding generating function |Blj . Thus, one may hope that extending 
Kuperberg's work to the 8VS0S model would lead to new applications of statistical 
mechanics to combinatorics. That is precisely the object of the present study. 

Apparently, the first step in this program is to generalize the Izergin-Korepin 
formula to the 8VS0S model. In a recent paper [R] we found such a generalization. 
In the trigonometric limit (which is intermediate between the six-vertex model 
and the general elliptic 8VS0S model), we used it to obtain a closed formula for a 
special case of the three-colour partition function, see f l2.5p . In the present paper, 
we consider the general case. Although there seems to be no very simple formula 
for the general three-colour partition function, we can express it in terms of certain 
special polynomials, which have remarkable properties and deserve further study. 
To obtain these results has not been straight-forward; in particular, we have not 
been able to work directly with the explicit formulas from [R] . Rather, we combine 
a simple consequence of those formulas with several further ideas. 

The plan of the paper is as follows. In ^ we describe the three-colour model 
with domain wall boundary conditions, and its relation to ice graphs and alternat- 
ing sign matrices. We refer to the states of the model as three-coloured chessboards. 
In ^ we state our main results in elementary form. ^ contains preliminaries on 
theta functions and the 8VS0S model. In ^ we prove our first main result. 
Theorem 13. ![ which expresses the three-colour partition function Z^*^ in terms of 
special polynomials g„ and r^. This is a rather easy consequence of results in |^. 

We then turn towards an alternative way of expressing Z^'-^ . In ^ we intro- 
duce a function $„, which provides a one-parameter extension of Z^'-' . For the 
six-vertex model, a similar function appears in Zeilberger's proof of the refined 
alternating sign matrix conjecture |Z2j . In §3 is generalized to a multivariable 
theta function \E'„. These functions play a similar role as Schur polynomials do 
in Stroganov's proof of the alternating sign matrix theorem [St] (see also [D]). 
However, while the Schur polynomials are instances of the six-vertex partition 
function (with the crossing parameter a cubic root of unity), the function \E'„ is 
not directly related to the 8VS0S partition function. The function \E'„ has two 
important properties, the first being a determinant formula reminiscent of the 
Izergin-Korepin determinant. The second property is a symmetry with respect to 
inversion of each variable, meaning that it naturally lives on the Riemann sphere 
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rather than the torus. As a consequence, can be identified with a symmetric 
polynomial Sn- In §Hl we specialize the variables in 5*^, obtaining certain two- 
variable polynomials P„ and one- variable polynomials p„. Using minor relations 
for the determinant defining \E'„, we obtain recursions for these polynomials, which 
can be used to derive many further properties. In ^ we return to the three-colour 
model, expressing Z^'-^ in terms of the polynomials pn- This result has combina- 
torial consequences. For instance, for three-coloured chessboards of fixed size, we 
can compute the maximal and minimal possible number of squares of each colour, 
see Corollary 13.31 

Finally, in §101 we study the thermodynamic limit n — )■ oo. Using non-rigorous 
arguments, we are led to an explicit formula for the free energy of the three- 
colour model with domain wall boundary conditions (Conjecture 13.141) . From the 
viewpoint of physics, this is the main result of the paper. To prove Conjecture l3.14l 
rigorously is an interesting problem, which we expect to be rather difficult. For 
the six-vertex model with domain wall boundary conditions, a rigorous analysis 
has been done only recently [BFl IBLTI [BL2l |BL3] . 

Note added in proof: Immediately after seeing an earlier version of the 
present paper, Vladimir Bazhanov and Vladimir Mangazeev sent me an interesting 
conjectured recursion for the polynomials p„. They have also obtained a similar 
conjecture for P^. These recursions give a much faster way of computing the 
three-colour partition function than those obtained in the present paper. Both 
conjectures can be found in jBM4j . where the authors stress the resemblance to 
polynomials occurring in their analysis of the eight-vertex model |BMlt IBM21 
IBM3] . It would be interesting to investigate if that relation can be made precise; 
for instance, that may give a link between the present work and the Painleve VI 
equation. 

Acknowledgements: This work was partly carried out while participating 
in the programme Discrete Integrable Systems at the Isaac Newton Institute for 
Mathematical Sciences, and I thank the institute and the programme organizers for 
their support. I would also like to thank Vladimir Bazhanov, Vladimir Mangazeev, 
Jacques Perk and, in particular, Don Zagier for their interest and for valuable 
comments and suggestions. Finally, I thank the anonymous referees for useful 
comments. 

2. Three-coloured chessboards 

We will refer to a state of the three-colour model with domain wall boundary 
conditions as a three-coloured chessboard. Fixing n, consider a chessboard of size 
(n-|-l) X (n + 1). The squares will be labelled with three colours, which we identify 
with the three residue classes 0, 1,2 mod 3. We impose the following two rules. 
First, vertically or horizontally adjacent squares have distinct colours. Second, 
the north-west and south-east squares are labelled and, as one proceeds away 
from these squares along the boundary, the colours increase with respect to the 
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cyclic order 

0<1<2<0. (2.1) 

In particular, the north-east and south-west square are labelled n mod 3. As 
an example, when n = 3 there are seven three-coloured chessboards; these are 
displayed in Figure [H 



Figure 1. The seven three-coloured chessboards of size n = 3, where is 
pictured as black, 1 as red and 2 as yellow. 




We will briefly explain the bijections to alternating sign matrices and ice graphs 
mentioned in the introduction. Let ( ^ d) be a 2 x 2-block of adjacent squares from 
a three-coloured chessboard, and choose representatives of the residue classes so 
that adjacent labels differ by exactly 1. For instance, if the original block is (21); 
we may choose the representatives ( 1 1 )• Having made such a choice, we contract 
each block to the number {h-\- c — a — d) /2, obtaining an n x n- matrix with entries 
— 1, 0, 1. For instance, from the last chessboard in f??) we obtain 




This gives a bijection from three-coloured chessboards to alternating sign matrices. 

To obtain the bijection to ice graphs, we draw an arrow between any two adja- 
cent squares in such a way that the larger label, with respect to the order (12. ip . 
is to the right. For instance, the last chessboard in (??) corresponds to the arrow 
configuration 



(2.2) 
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The result is a directed graph, where each internal vertex has two incoming and 
two outgoing edges. Considering vertices as oxygen atoms and incoming edges as 
hydrogen bonds, this can be viewed as a model for a two-dimensional sheet of ice. 
We are interested in the generating function 

^colour 

chessboards squares 

^ 2.3 

fco+fci+fc2 = (n+l)2 

where N{ko, ki, ^2) denotes the number of three-coloured chessboards with exactly 
ki squares of colour i. In physics terminology, Z^'" is the partition function of the 
three-colour model with domain wall boundary conditions. 

It is difficult to study Z^^ by direct methods. Indeed, to compute ^^""(1, 1, 1) = 
An was an unsolved problem for more than a decade. In [R], we generalized that 
enumeration using the trigonometric 8VS0S model. Namely, we found a closed 
expression for ^^^(^0,^1,^2) when 

This surface can be parametrized by tj(A, yu) = — Acu*)^, where, as throughout 
the paper. 

By homogeneity, we may take yU = 1. Then [Rl Cor. 8.4], 
1 1 1 



3C 



(1-A)3' (1-Au;)3' (1-Aw2)3^ 

= ^'"(^-A3)l"ir'^' (-4-(l + -"^^) + (-D-^nC-A), (2.5) 
where is as in (11. ip and 

C„^2-5...(3n-l) 
" l-4---(3n-2) ^ ^ 

is the number of cyclically symmetric plane partitions in a cube of size n jX] . 
Although (12. 4 p is a strong restriction, it is sufficient for computing the moments 

^ N{ko,ki,k2)k^, 2 = 0,1,2, 

feo+fci+fc2 = (n+l)2 



see lEl Cor. 8.5]. 
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3. Statement of results 

3.1. Polynomials g„ and r„. In this Section, we state our main results. We 
begin with the following fact. As we will see in ^ it is a rather straight-forward 
consequence of results in jR] . 

Theorem 3.1. Let 

T=T(to,t,,t,)= ^^°^^+;f (3.1) 

Then, there exist polynomials qn andvn such that, forn = mod 3, Z^*" (to, ^i, 1^2) 
equals 

while for n = 1 mod 3 it equals 

(-1)"^^(W,)^ ( ^g.(T) -2-("°'^<^)^^^^l±^^^l±^,4T)) 

V ^^2 ^2 / 

and for n = 2 mod 3 it equals 

(-1)"+^(W2)"^ L(T) -2x("°-)^2!l±|!!i±^,„(T)) . 

\ tot2 J 

Here, x(true) = 1, x(false) = 0. 

The first few instances of the polynomials g„ and r„ are given in Table [H 
Note that, in each case, Z'^ is symmetric in the two variables t_„±i mod 3- This 
is explained by the fact that reflection in the vertical (say) axis, followed by in- 
terchanging the colours — n ± 1 mod 3, defines an involution on three-coloured 
chessboards. Theorem 13.11 shows that Z'^ is very nearly symmetric in all three 
variables, being a linear combination of two symmetric polynomials, where the 
coefficients are polynomials in t_„ mod 3 of low order. Moreover, the symmetric 
polynomials depend only on the second and third elementary symmetric poly- 
nomial, being independent of to + + '^2- 

Table 1. The polynomials qn and r„. 



n g„(x) r„(a;) 









1 


1 


1 





2 


1 


1 


3 


1 


1 


4 


X + 3 


X — 3 


5 


— 4x + 6 


X + 6 


6 


x^ - 2x + 40 


x^ - Sx^ + 20 


7 


- lOx^ + ISx^ H 


V lOOx + 50 x^ + 75x - 50 
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The following result seems much deeper than Theorem I3.1[ We need consider- 
able preparation for its proof, which is given in §9.21 

Theorem 3.2. The polynomials and r„ are monic. Moreover, their degrees are 
given by 

,2 



deg(g„) + 1 = deg(r„) 
deg(g„) = deg(r„) + 1 
deg(gn) = deg(r„) 
deg(g„) = deg(r„) 



n 
12' 



12 ' 

2 __4 



n 



12 ' 
2 9 



n 



12 



n = mod 6, 
n = ±1 mod 6, 
n = ±2 mod 6, 
n = 3 mod 6. 



As an application, we can determine the maximal and minimal number of 
squares of each colour. These bounds restrict the counting function N introduced 
in (12.31) to an equilateral triangle. We can also explicitly evaluate the restriction of 
N to the boundary. To formulate the result, we introduce some notation. Fixing 
n, let 

{N{x,y,z), n = mod 3, 
N{y,z,x), n = l mod 3, 
N{z,x,y), n = 2 mod 3, 

so that N{x, y, z) = N{x, z, y). Moreover, let 

'(5n2 + 8n)/12, n = 0, 2 mod 6, 

(5n2 + 8ra + 11)/12, = 1 mod 6, 

(5^2 + 8n + 3)/12, n = 3, 5 mod 6, 

^ (Sn^ + Bra + 8) /12, = 4 mod 6, 

'(^2 + 4n)/6, n = 0, 2 mod 6, 

{n^ + 4n + 7)/6, n = 1 mod 6, 

{n^ + 4n + 3)/6, n = 3, 5 mod 6, 

^(?T,2 + 4n + 4)/6, n = 4 mod 6, 



M 



5n^ + Bra + 11 
12 



m 





'I, 


n = 0, 2 mod 6, 




-2, 


n = 1 mod 6, 


= < 


0, 


n = 3, 5 mod 6, 






n = 4 mod 6, 




H 


(n2-l)/4, n odd. 


_ 4 




n2/4, n even 



(3.2) 



8 HJALMAR ROSENGREN 

Note that 2M + m + e = {n + if. Let 

A = {{x,y, z) e Z^; X + y + z = {n + if, x <M + e, y, z < M}. 

Then, A is an equilateral triangle, with S lattice points on each side. We denote 
its corners by 

P = {m + €,M,M), Q = {M + €,m,M), R={M + e,M,m). 

Corollary 3.3. The convex hull of the support of N is equal to A when n is odd 
and A \ {P} when n is even. In particular, the maximal number of squares of 
each of the colours —n ± 1 mod 3 is equal to M , and the minimal number of such 
squares is m. The maximal number of squares of colour —n mod 3 is M + e, and 
the minimal number is m + e if n is odd and m + e + 1 if n is even. Moreover, the 
restriction of N to dA is given by 



N 



kP + {6- k)Q 



N 



kP + {6- k)R 
~5 



5-1 
k 

5 
k 



n even, 



n odd, 



N 



kQ + {5- k)R 



6 
k 

6-2 
k 



+ 



6-2 
k-2 



n even. 



n odd. 



where < k < 6 . 



It is straight-forward to derive Corollary 13.31 from Theorem [321 some details are 
given in §9.21 

As we explain at the end of ^ it follows easily from Theorem 13.11 that the 
polynomials g„, r2„ and 2r2„+i have integer coefficients. However, the following 
fact is not obvious, see §9.3[ 



Proposition 3.4. The polynomial r2„+i has integer coefficients. 

This result has a simple combinatorial meaning. If n = 3 or 5 mod 6, Theo- 
rem 13.11 expresses Z^*" as a sum of a symmetric polynomial and a polynomial with 
even coefficients. Thus, the function 

{ko,ki,k2) ^ N{ko,ki,k2) mod 2 

is symmetric. Similarly, when n = 1 mod 6, 

{ko, ki, k2) ^ N{kQ, ki, ^2 - 2) mod 2 

is symmetric. In the notation of Corollary 13. 3[ these facts can be stated as follows. 

Corollary 3.5. When n is odd, N mod 2 is invariant under the action of S3 as 
the symmetry group of A. 
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To illustrate Corollaries 13.31 and 13.51 we give two examples. When n = 4, the 
non-zero values of are 



(Ar(6 + z,6 + j,13-z-j))J,.=o 



/ 1\ 

4 3 
6 6 3 
4 6 1 
yi 3 3 1 

Since n is even, the lower right corner is missing from the support of A^, and all 
boundary entries are binomial coefficients. When n = 5, the non-zero values are 



(iV(8 + ^,20-2-j,8+j)) 



i,j=0 













1\ 










4 


6 








7 


18 


15 






8 


12 


36 


20 




7 


12 


36 


40 


15 


4 


18 


36 


40 


24 


6 


\1 6 


15 


20 


15 


6 


1/ 



In this case, the diagonal entries are numbers of the form + (^Ig)- Note the 
symmetric distribution of the odd entries, which is pecuhar to the case of odd n. 

3.2. Polynomials p„ and The polynomials g„ and r„ are closely related to 
a third class of polynomials, which we denote pn- In the following result, which is 
proved in §9.2[ we use the notation 

p(x) = x'^'^Ppil/x); (3.3) 

that is, p denotes the polynomial obtained from p by reversing the coefficients. 

Theorem 3.6. There exist polynomials pn of degree n{n + l)/2 such that, for n 
odd, 



p„.i(C) - C'^'pn-ilC) = (1 - C)(C' + 4C + 1)^ 



C(C + i)' 



Pn-l(C) - C^Pn-l(C) = (1 - 0(1 + C)'(C' + 4C + l)^r„ 

whereas for n even, 

Pn-i(C) - C"*'p„_i(C) = (1 - C)(C' + 4C + l)^r, 



2(C2 + 4C + 1)3 

c(c+ir 



2(C2 + 4C + 1)3 

c(c+ir 



P._i(C) - C^Pn-i(C) = (1 - C')(C' + 4C + i)"^gn 



2(e + 4C + l)3 

c(c+ir 

2(e + 4C + l)= 
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Table 2. The polynomials pi 



-1 1 

1 

1 3C + 1 

2 5C' + 15C' + 7C + 1 

3 |(35C^ + 231C^ + 504C^ + 398C^ + 147(2 + 27C + 2) 

4 i(63C^° + 798C^ + 4122C^ + 11052C^ + 16310C^ + 13464C^ + 6636C^ 

+2036C^ + 387(2 ^ + 2) 

5 i(231C^^ + 4554C^^ + 39468C^^ + 196922(^2 + 622677C^^ + 1298446C^° 

+1816006C^ + 1726302C® + 1140593C^ + 535478C^ + 181104C^ 
+44134C^ + 7603C^ + 882(2 ^_ g2( + 2) 

6 1(1716(2^ + 50193(2° + 673530C^^ + 5484050C^*^ + 30199260C^^ 

+118703208C^^ + 342834244C^^ + 738954900C^^ + 1198556100C^^ 
+1470762970(^2 + 1373623128C" + 984509064C^° + 546520100C^ 
+236837400C^ + 80476380C^ + 21422188C^ + 4430904C^ + 699405C^ 
+81550C^ + 6630(2 ^ ggg^ ^ 



The shift in n is introduced for convenience. See Table |2] for the first few 
instances of the polynomials Pn- 

To indicate the meaning of Theorem 13.61 we solve for pn, obtaining for n even 

.Cir.0He.^C.r)^,,..{,^^ii^,)) (3.4) 

and for n odd 

A priori, the right-hand sides are polynomials of degree n{n + 2)/2, {n + 1)2/2, 
respectively. The degree bound n{n + l)/2 imposes relations between high coeffi- 
cients of the polynomials g„ and r„. In the notation of Corollary 13. 3[ this implies 
relations between values of N close to the boundary dA; the details will not be 
worked out here. 

An important property of the polynomials pn is that they appear as solutions 
to certain linear equations. 
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Theorem 3.7. Consider the polynomial equation 

AX-BY = C, (3.5) 

where 

A = (c + i)V-i(c), B = c^'p^{c), 

/ A\ 1+X(" odd) 



n(n+l)/2 



c = (i + 2C)^+^^"^^^°^(^i + ^j Pnicy. 

Then, this equation has a solution {X, Y) such that X = Pn+i and degY = degA. 

Theorem 13.71 follows from Proposition 18.31 and Corollary 18.101 Although we 
have not been able to prove it, we believe that the polynomials A and B are 
always relatively prime. In fact, we believe that p„ is irreducible over Q for 
n > 2. Assuming that this is the case, (13.51) has a unique solution (Xq, Yq) with 
degYo < degA. The solutions with degF < degA then have the form {X,Y) = 
{Xq + XB, Yq + XA), A G C. We can specify the solution X = Pn+i in this space 
by declaring that the leading term is 

2-m ^2^+2' 

see Proposition 18.91 Then, Theorem 13.71 gives a recursive procedure for construct- 
ing the polynomials p„. This gives a comparatively fast method for computing 
the partition function Z^'-' . (As was mentioned at the end of the Introduction, a 
much faster method has been suggested by Bazhanov and Mangazeev |BM4j .) 

The polynomials p„ have coefficients in see Corollary l9.8l for a more precise 

statement. Moreover, the following facts seem to be true. 

Conjecture 3.8. The polynomials Pn have positive coefficients. 

Conjecture 3.9. The polynomials pn are unimodal in the sense that, if 

n(n+l)/2 

PniO= Yl ""^^^^ 
k=0 

then, for some N, 

oo <«!<•■■ < dN-i < aN > aN+i > ■ ■ ■ > a„,(„+i)/2. 

Using Theorem 13. 7t we have verified Conjectures 13.81 and 13.91 up to n = 16. For 
these values, the maximal coefficient occurs at 

{n(n + 2)/4, n even, n < 16 

{n + 1)2/4, n odd, n<7 

(n- l)(n + 3)/4, n odd, 9 < n < 15. 

Numerical experiments suggest that the zeroes of pn form quite remarkable 
patterns, see Figure [2l 

We have some partial results on the real zeroes, see §8.41 for proofs. 
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Figure 2. The 105 complex zeroes oi pu. 



-1 




--2 

O 

o --3 

Conjecture 3.10. The polynomial pn has exactly [{n + l)/2] real zeroes. 

Conjecture 3.11. The polynomial pn does not vanish in the interval —2<(< 
-1/2. 

Proposition 3.12. Assume Conjecture I3.10[ Then, all real zeroes of pn{C) are 
simple, and contained in the interval — 1/2 <(< if n is odd and in ( < —2 if n 
is even. Moreover, the real zeroes of p2n interlace the real zeroes of p2n+i, which 
in turn alternate left of the real zeroes of p2n+2- In particular. Conjecture 13.101 
implies Conjecture 13 .111 

The polynomials Pn can be embedded in a more general family of polynomials 
Pn{x, C), which are of degree n in x and degree [n{n + 2)/2] in To be precise, 

P„(l + 2C,C) = (l + 2C)Wpn(C)- 

We give the first few instances of Pn in Table [3l 

The polynomials Pn satisfy a three-term recursion of the form 

A-nPn+l = BnPn + CnPn-li 

where An = An{x,() and C„ = C„(x, C) are quadratic polynomials in x, while 
Bn = Bn{x,() is cubic in x, see Proposition 18.61 This is reminiscent of the 
recursion satisfied by orthogonal polynomials (where An and C„ are constants 
and Bn is linear). It also seems that P„ resemble orthogonal polynomials with 
respect to their zeroes. Indeed, the following fact is proved in §8.4[ 

Proposition 3.13. Assume Conjecture 13.111 (or, in view of Proposition I3.12| 
Conjecture 13.10)) . Then, if —2 < ( < —1/2 and ( ^ —\, all zeroes of Pn{xX) are 
simple and positive. Moreover, the zeroes of Pn+i interlace those of Pn. 
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Table 3. The polynomials P„. 

n PnjxX) 

1 

1 x + C 

2 l((C + 2)(3C + l)x' + C(C + 3)(3C + l)x + eiC + 3)(2C + 1)) 

3 i^(C + 2)(5C=^ + + 7C + l)x'^ + C(C + 4)(5C^ + ISC' + 7C + l)^;' 

+c'(4c + i)(c' + re + 15C + 5)x + C''(2c + i)(c^' + 7e + i5c + 5)) 

4 I ((C + 2)'(35C^ + 231C^ + 504C^ + 398C^ + 147C' + 27C + 2)x^ 

+C(C + 5)(C + 2)(35C^ + 231C^ + 504C^ + 398C^ + 147C' + 27C + 2)x^ 
+3C'(10C^ + 139C^ + 790C'^ + 2245C^ 

+3232C^ + 2245C^ + 790C^ + 139C + lO)^^ 
+C^(2C + 1)(5C + 1)(2C^ + 27C^ + 147C^ + 398C^ + 504(2 ^ 231C + S5)x 

+C^(2C + 1)2(2C^ + 27C^ + 147C^ + 398C'^ + 504(2 ^ 231C + 35)] 

Computer calculations suggest that, if —2 <(<—!, the zeroes are in fact 
contained in the interval x > 1, while if —1 < C < they are contained in 

< a; < 1. 

3.3. Symmetric polynomials S'„. The key to the proof of most our results 
is that the polynomials P„ and Pn can be obtained as specializations of certain 
symmetric polynomials Sn of 2n + 1 variables. To define them, we introduce the 
elementary polynomials 

/(x) = (C + 2)x-C, g{x)^{x-l){x-C), h{x)^x'{x-{2C + l)), (3.6) 
and then let 

1 (fix) m f(z)\ 

F(x,y,z) = det q(x) q(y) qiz) 

(y - - x)iz - y) [II] ill] (3.7) 

= (C + 2)xyz - C{xy + yz + xz + x + y + z) + C(2C + 1), 

= (C + 2)xy{x + y)- ({x^ + y') - 2{e + K + l)xy + C(2C + l){x + y). (3.8) 
In this notation, 

Sn{Xli • • • 1 Xfi-i Vl-i • • • 1 Vm -^) 

_ ^^.,Gfe.*) f^^filWlfiy (3.9) 
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the dependence on ( being suppressed from the notation. Though it is not appar- 
ent from this definition, Sn is a symmetric polynomial of all 2n + 1 variables, see 
§7.2[ Moreover, 

5„(.,2C + 1.....2C + l,^.....^) 

n 

r (n-l)^ 1 _ 2 

= (-i)("iV(i+cr'(i + 2C)[^](i + 0^ ' ' p„{x, 

see Proposition 18.11 

We mention that the special case ( = —2 of S'„ is an orthogonal character: 



0, 



Sn{Xi, . . . ,X2n+l) "2" xln,n,n-l,n~-l,...,l,l,0)('^'^^ ' ' ' '^2n+l), 

C— 2 

where Xj = —(l + ti + t^^); see Theorem 17.171 for a more general result. For 
general (, Sn can probably be interpreted as an affine Lie algebra character; see 
Remark 17.41 

3.4. Thermodynamic limit. All results mentioned so far concern properties of 
the partition function for fixed n. However, from the viewpoint of statistical 
mechanics, the main problem is to investigate the asymptotics as n — )■ oo. A 
particularly important quantity is the free energy per volume, which we identify 
with the limit 

f{to,'ti,t2) = lim . (3.10) 

n— )-oo n 

We propose the following explicit expression for /. In ^ ITO] we present a formal 
derivation of this result; to give a rigorous proof is presumably quite difficult. 

Conjecture 3.14. Assume that the parameters ti are all positive, and let ( be any 
positive solution of 

_ (tpti + tot2 + t,t2f _ + 4C + 1)' , . 

(W2)2 c(c + i)^ ■ ^ ■ ^ 

Then, the free energy per volume is given by 

f{toM,t2) = ^ log(totit2) + log [ (C + 2)i(2C + l)i \ _ 
3 y 23Ct2(C + 1)3 ) 

Note that, by the arithmetic-geometric inequality, T > 27. Using that 
2(C^ + 4C + 1)=^ (C-1)^(C + 2)(2C + 1) 

it is easy to check that (13. lip always has two positive solutions, except in the case 
T = 27, when C = 1- Moreover, the two solutions are related hj ( ^ which 
does not change the right-hand side of ( 13.12^ . 
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As an example, when to = ti = t2 = C = f l3.12p gives 

'3^3' 
T 

which agrees with the known asymptotics |BF] 



/(1, 1,1) = log 



Zf(l,l,l) = A„~Cn-^ l^^j , n^oo. (3.14) 

One should compare Conj ecture 13 . 1 41 with Baxter's result for periodic boundary 
conditions |B1] . In that case, one still has 

f{to,tiM) = ^log(toM2) + logVr, 

with W a function only of T. Baxter gives the formula 

64(l-9t^)§ 
27(l + t)3(l-3t)' 

where t is the unique solution of 

T il- 3^2)3 



27 1 - 9^2 

such that 0<)f:<l/3. It is straight-forward to check that t = (C — 1)/3(C+1), 
where C, is the unique solution of (I3.1ip such that C ^ 1- This gives 

2ic^(C + i)t 



w = w. 



(2C + i)i ' 

which is manifestly different from our conjectured formula 

(C + 2)f(2C + l)t 



DWBC 



2fC^(C + l)i 

for domain wall boundary conditions. Note the intriguing relation 

2 

W^DWBC(C) 



4. Preliminaries 

4.1. Theta functions. We will work on the multiplicative torus / {z = pz}, 
where C* = C \ {0} and < |p| < 1. We introduce the theta function 

00 

e{x;p) = Y[{l-p^x){l-p^^^/x). 

j=0 

We often employ condensed notation such as 

9{ai, . . . ,a„;p) = 6'(ai;p) ■ ■ ■6'(a„;p), 
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e{xy^;p) = e{xy;p)e{xy ^■,p). 
The most fundamental relations for the theta function are 
9{px;p) = 9{x~^;p) = — x~^6'(x;p), 
together with the addition formula 

9{xz^ ^yw^] p) — 9{xw'^ , yz^]p) = — 9(xy^, zw^;p). (4.1) 

z 

Other elementary identities that we will use include 

9{x'^;p'^) = 9{±x;p) = 9{±x, ±px;p^), 
from which one can derive 

^(-l,±p;p2) = 2, (4.2a) 
9{-uj, ±puj; p^) = -Lo^, (4.2b) 

where u = e^'^*/^. 

We will use the following terminology from |RSc] . 

Definition 4.1. A holomorphic function on C* is called an An-i theta function 
of nome p and norm t if it satisfies 

It is called a BCn theta function of nome p if 

Finally, it is called a Dn theta function of nome p if 

Lemma 4.2 ( [RScj Lemma 3.2]). A function f is an An-i theta function of nome 
p and norm t if and only if it can he factored as 

f{x) = C9{aix, . . .,anx;p), 

where C G C and ai G C* with ai ■ ■ ■ a„ = t. It is a BCn theta function of nome 
p if and only if 

f{x) = C9{x, ±y/px; p)9{aix^ , . . . , an-ix^; p) , 

where C E C and G C* . Finally, it is a Dn theta function of nome p if and only 
if 

f{x) = C9{aix^, an-ix^]p), 
where C G C and G C*. 
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The terminology is motivated by Macdonald's theory of affine root systems |Mlj . 
In each case, the Macdonald identity for the affine root system R is equivalent to 
evaluating a determinant det(/i(xj)), where fi runs through a basis in the space 
of R theta functions, see |RSct Proposition 6.1]. We will need two special cases, 
which are both classical theta function identities. For R = Ai, the Macdonald 
identity can be written 

9(ax,bx;p) = — (9(—pa/b,—abx'^;p'^) — bx9(—a/b, —pabx'^;p'^)) ; (4.3) 

"[P',P) 

this can also be obtained from ( 14.11) . The Macdonald identity for BCi is Watson's 
quintuple product [W] 



where 



e{x, ± Vpx; p) = {e{-px^; p3) - xei-px-""; p"^)) , (4.4) 

[Pi P)oo 



{p;p)oo = llii-p^). (4.5) 

i=i 

Finally, we recall some classical facts on uniformization. The Riemann surface 
S = C* /{z = pz = z~^} is a sphere, the analytic automorphisms ^ : S —¥ CU{oo} 
being given by 

^^''-^eibz^;p)' 

where C G C* and a, b are distinct in C*/ {z = pz = z~^}. Accordingly, any Dn+i 
theta function / can be written 

f{z) = d{bz^;prP{az)), 

where P is a polynomial of degree at most n. We will refer to the passage from 
/ to P as uniformization. When / is factored as in Lemma 14. 2[ the explicit 
uniformization is 

^± ^ ^±.^^ _ ai---anO{baf,...,bat;p) 
ti[a,z ,...,a^z ,p)- c-a-e{ba^-pY 

X e{bz^;pr{a^) - e(ai)) ■ ■ ■ (e(^) - ^m). (4.6) 

This follows immediately from the case n = 1, which is equivalent to f l4.ip . 

4.2. The 8VSOS model. The 8VS0S model is an ice model, so for domain 
wall boundary conditions states can be identified with three-coloured chessboards. 
Fixing the size of the chessboards to (n + 1) x (n + 1), the model depends on 
parameters 

a;i,...,x„,?/i,...,?/„,A,pe C*, \p\ < 1. 

In general, there is also a crossing parameter g, but for our purposes it can be 
fixed to a cubic root of unity. 
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We will assign a weight to each 2 x 2-block of adjacent squares. These blocks 
can be viewed as entries of an n x matrix, and are given coordinates 1 < i, j < n 
in a standard way. A block (° ^) with coordinates is then given the weight 

Ra~_ltai>'^''^^^/yj)^ where 

Here, ± is a short-hand for ±1 mod 3. 
The partition function is now defined as 

Zl^^o^{xi,...,Xn;yu...,yn;Kp)= ^ JJ weight (block). 

chessboards blocks 

This differs slightly from the normalization used in [R]. We have 
Zf soS(^. y. x^p^ = {y,... y^y^Z^ix- y- A), 

where Z„ is as in |Rl §7]. 

In the specialization Xi = u, yi = 1, the 8VS0S model reduces to the three- 
colour model with parameters 

More precisely, in the case of domain wall boundary conditions, 
111 



0(Aa;-p)0(A3;p3).^+2n+2^n (c^, • • • , 1, . . . , 1, A, p) . (4.8) 

See [HI §8] for the trigonometric case p = 0; the discussion there carries over 
verbatim to general p. 

We need to mention the recursion 



Z^^^^{xi, . . . ,Xn;yi,. . .,yn; A,p) 

;p) 111=2 ^jyi^'^/Vk, 
^(Aw"-i;p)^(w;p)2'^-2 



Luxi=yi 



<^ 777\ — ~i TTTT \2r,-'2 [^2, ■ ■ ■ , ^n, y2, ■ ■ ■ , yn, ^, P) , 



(4.9) 
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see e.g. [R| Lemma 3.3], and the crossing symmetry 
^8VSOS(^2/^^^ . . . , a;Vx„; IM, . . . , l/y„; u?^ I\ p) 



x;y;X,p), (4.10) 



^(Aw";p)X" 
where 

X = Xi---Xn, Y = yi---yn. 

The equation fl4.1Qp can be derived from a corresponding symmetry of the Boltz- 
mann weights, and is also apparent from the exphcit formulas for the partition 
function given in [R]. 

The main result of [R] is the explicit expression 



X 



e{u; pY e{r. p)"r"+i^(XA7w"/r ; p) 

n"j=i yp{xi/yj,ujXi/yj;p) 
ni<i<i<n Xjyje{xi/xj,yi/yj;p) 

Sc^„^n} 0{\u'^-\S\-p) l<r,3<n\e{xl/yf,p) 



(4.11) 



where 



a;iCi;, i E S, 
Xi, i ^ S 



and 7 is arbitrary. This can be viewed as an analogue of the Izergin-Korepin 
formula for the six-vertex model. We will not work with this formula directly, 
though we need the following immediate consequence. 

Proposition 4.3 ((Rj Corollary 5.4]). As a function of X, 

e{Xu-^\Xu-^';p)Zl'''^'{x-y-X,p) 

is analytic on C* . More precisely, it is an Ai theta function of nome p and norm 

By dSHD, it follows that 



^(Aw2^Aw"+i;p)2 " \e{X;py^' e{Xuj;pY' e{Xio^;py 

is analytic in A G C*. This remarkable fact is a key result for the present work. 
We do not know how to prove it except as a consequence of (14. lip . 

We also need the following functional equation, which was recently obtained by 
Razumov and Stroganov. 
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Proposition 4.4 ([RS]). Let 

Fn{x; y; A, p) = e{Xoo^^\ Xoo^^'; p)A{x, y; ^^^(cux; y; A; p), 

where 

A{xi, . . . ,xn;p) = Y\. Xj9{xi/xj;p). (4.13) 

Then, 

2 

^Fn{uj''xi,X2, . . .,Xn;y;uj~''X,p) = 0. 

fc=0 

Consider F„ as a function of A. By Proposition 14.3^ it is an Ai theta function 
of nome p and norm u'^Y/X. By (14.31) . it can be decomposed as 

F^{x;y;X,p) = Fi'\x;y;p)e{-io^X'Y/X;p') + Fi'\x;y;p)Xe{-pio^X'Y/X;p'). 

The following result is then immediate from Proposition 14. 4[ 

Corollary 4.5. For i = 0,1, 

2 

^F»(w'=Xi,X2,...,x„;y) = 0. 

fc=0 

5. Proof of Theorem 13.11 

Consider (I4.12p as a function of A. By Proposition 14.31 it is an Ai theta function 
of nome p and norm w". Since the space of such functions is spanned by ^(Ao;^"; p)^ 
and 6{Xijj'^"'~^^, Acj^""*"^; p), there exist functions X„ and y„ such that 

^(A3;p3)n^+2n+3 



eiXu^^Xu'^+^-pf 

= X„(p)^(Au;2";p)2 + F„(p)^(Au;2"+\ Xu^'^+^-p), (5.1) 

where ti are as in (14. 7p . We will see that this decomposition corresponds to 
Theorem 13.11 the p-dependence being encoded in the polynomials g„ and r„. 
We first give some preliminary results. 

Lemma 5.1. There exists a function p ^ t{p) such that, under the parametriza- 
tion (gZD, 

i + i + l = r(p)^(A3;p3). (5.2) 

to tl 12 

Moreover, the quantity T (to, ti,t2) , defined in (13. ip . equals rijpY . 

Proof. Each of the functions 1/ti = 6{Xu^;p)^, as well as ^^(A^;p'^), are A2 theta 
functions of nome p and norm 1. Thus, the first statement can be reduced to the 
trivial verification that ^jl/^i vanishes at the points A = l,u,u'^. The second 
statement follows using 

_ 1 

^otit2 - J^^^s- (5.3) 
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□ 



It was pointed out to us by Don Zagier that 



rip) = 3(1 + 9p^^^^ ) , (5.4) 



where we use the notation (14. 5p . We also have (see the remark after Lemma [9. ip 

r{pf = 27(1 + 27p ^f'^2f ) ■ (5-5) 



This is consistent with (15. 4p in view of the identity 



{p;p)lo J {p;py^ ' 

which can be found in Ramanujan's notebooks |Be[ p. 345]. We mention that the 
function 

ip^;p^y^ 
(p; pY^ 

is well-known in the theory of modular forms. It is automorphic under the group 

ro(3) = I rf) ^ ^ = 3 

(acting by r (ar + 6)/(cr + d), where p = e^'^*'^) and in fact generates the field 
of all such functions [SI Thm. 21]. 

Lemma 5.2. Let f be a homogeneous rational function in three variables. Suppose 
that, under the parametrization ( 14. 7p . /(to, ^1,^2) = 0. Then, / = 0. 

Proof. It is enough to show that the Jacobian of the map (A, /x, p) 1— )■ (yuto, /^^i, ^^^2) 
does not vanish identically. Clearing denominators, the Jacobian is proportional 
to 

/ e{x-p) 9',{x-p) ff^{\-p) 

det e{\u-p) uj9'^{\lo;p) e'p{\u;p) 
\e{Xu';p) u'e'^iXu'^p) e'^{\u';p)^ 
Using that 

e{x; p) = (1 - x) (1 - (x + x'^)p + (1 - a; - x'^)p^) + 0{p^), 
one can check that the determinant is 

3V3.(H-2A3)(l-A^^)^^ 

□ 

Lemma 5.3. Let f be a Laurent polynomial in three variables, homogeneous of 
degree 0. Suppose that, under the parametrization (14.71) . fitQ,ti,t2) is independent 
of A. Then, f is a polynomial in T. 
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Proof. By the change of variables A uX, it is seen that /(to, ti, h) = f(ti, t2, to)- 
A priori, this holds when ti are as in (14. 7p . but by Lemma \572\ it is valid in general. 
Similarly, A i— A~^ gives /(to, ti, ^2) = /(to, t2, ti). Since these two transformations 
generate S'3, / is symmetric, and can thus be expressed as 

/(to, tl, t2) = Cklmito + ti + t2)^(totl + tot2 + tit2)'(totit2)™'. 

fc+2/+3m=0, fe,i>0 

Introducing the function 

0(A,p) = (to + ti+t2)^^(A3;p3)3 

and using (15. 2 p and (15. 3 p gives 

A;>0 ^ ^ 2l+3m=-k, l>0 

We observe that, since 

the fcth term has a pole at A = 1 of order exactly k. In particular, if / is indepen- 
dent of A, then all terms with k ^ vanish. Thus, 

/(to,ti,t2)= C'oim(totl + tot2 + tit2)'(totit2)™, 

2l+3m=0,l>0 

which is indeed a polynomial in T. □ 

We are now ready to prove Theorem 13.11 Assume that n = mod 3. Then, 
(15. ip can be written 



(^0, tl, t2) - 7T7T^r-iV + 



^(A3;p3)"("+2) V^(A3;p3) ^(A;p)3 
(tot,t2)^ ( ''-{p)r{p)hUt2 ^ ^^^^^^^ 

\toti + tot2 + tit2 



where r(p) is as in Lemma [5.11 This implies 

V f \ f \ toti + tot2 + tit2 tiZ„(to, tl, t2) — to^„(ti, to, t2) fr n\ 

Xn{p)r{p) = 7 7 • (5.6) 

(totit2)^ + l ^1-^0 

Since the right-hand side vanishes when T = 0, it follows from Lemma [5.31 that 
Xn{p)T{p) is a polynomial in T divisible by T. Similarly, 

Y( \ _ ^ Zn(tQ,ti,t2) — Zn(ti,to,t2) 

(totlt2)^ ^1-^0 

SO Yn{p) is a polynomial in T. Writing 

Xn{p)r{p) = (-l)"+^Tg„(T), (5.7a) 
y„(p) = (-l)"2>^('^°<^'^V„(T), (5.7b) 
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we obtain Theorem 13.11 for n = mod 3. The cases n = ±1 mod 3 are similar 
and we do not give the details. For later use we note that in both cases 

X4p) = (-l)"+ig„(T), (5.7c) 

YnipMp)^ = (-l)"2>^("°'^^)Tr„(T). (5.7d) 
Finally, we comment on the claim made before stating Proposition 13. 4[ that g„, 
r„ and 2r2n+i have integer coefficients. Consider the case of g„ when n = mod 3; 
all other cases follow similarly. By (15. 6p and (I5.7ap . 

tlZnito, ti, t2) — tQZn(ti,to, ^2) 
ti — to 

V {totlt2j J 

The left-hand side is a symmetric polynomial in ti with integer coefficients. By 
|M2t 1.2.4], it can be expanded as an integer linear combination of elementary 
symmetric polynomials. By the identity above, the non-zero coefficients in that 
expansion are coefficients of the polynomial 

6. The function $„ 

In ^ we expressed the function (I4.12p in terms of the basis 
(^(Au;^p)^^(Au;2n+l^^^2n+2.^)^^ 

The first main idea for analyzing the partition function further is suggested by the 
trigonometric case p = 0. The relevant space is then the polynomials A -|- B\ -|- 
Cy? such that C = uj"'A. From (12. 5p . it appears that the most natural basis to 

use is not ((1 - XCO^'')^, (1 - Aw2"+l)(l - Aw2n+2)) ^^^^^^ (^^ ^ ^n^2^ 

elliptic analogue of the latter basis is {9{—ijj"'\'^;p'^),\9{—pu^X'^;p'^)) . As we will 
eventually see, this change of basis corresponds to expressing in terms of the 
polynomials Pn and Pn rather than g„ and r„. 

The next main idea is to consider a one-parameter extension of Z^'~^, which 
is given by Z^^^^ with A, p and Xi free but all other parameters fixed. This 
corresponds to incorporating Pn(C) the two-variable polynomial Pn{x, (). 

Combining these two ideas, we consider the function 

e{Xco-^\ Xu-^'- p)Z^^'''^{cot, co,...,oj;l,...,l;X,p). (6.1) 

By Proposition 14.31 and (14. 3p . it is a linear combination of 6{—u"'X'^/t;p'^) and 
X9{—poj^X'^/t;p'^), with coefficients independent of A. It will be convenient to 
write the corresponding decomposition as 

Zf^o^(a;t,a;,...,a;;l,...,l;A,p) ^ 



^(Aa;"+i,Au;'^+2;p)^(t;p)2"-i 
X (p"^t"^^(-a;"AVt;p2)l>„(t) -p-4^a;'"Ae(-pa;"AVt;p=')<l>„(t)' 



3n-4 9n-4 
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for n even, while for odd n we write 



78VSOS 



1 



3n — 3 9n — 5 



The functions and depend imphcitly on p but are independent of A. 

Lemma 6.1. The functions and are related by $„(t) = Moreover, 

1 



,6n-4-/-6n-3 



pi 



(6.2) 



Proof. This foUows easily from the fact that, as a function of t, (16. ip is an An-i 
theta function of nome p and norm w"/A, see [Rl Lemma 3.2]. □ 

As we will see, is uniquely determined by the following properties. 

Proposition 6.2. The function has the following properties: 

(i) $n is a BC^n-2 theta function of nome p^ ; 
2 

(a) ^a;'^$„(a;'^t) = 0, that is, $n(i) = f{t^) +tg{t^), with f and g analytic on 

k=0 

C*; 

(Hi) t = 1 and t = p are zeroes of ^n{t) of multiplicity at least 2n — 1; 
9{p;p^) 

Proof. We have already noted the quasi-periodicity (16.21) . For (i), it remains to 
show that $„(t"^) = -t"^$„(t). This is a special case of (I4.10p . 

Property (ii) is a special case of Corollary 14. 5 [ where we should note that 



(3n-2)/2 



(l-3n)/2 



A{x,y;p) 
Fi'\x-y) 



A{x,y;p) 



Xi=pt,X2- 



Xl=t, X2- 



= Xn=yi = 



=yn=i 



n even, 



n odd, 



where A is as in (14.130 and the hats indicate omission of Xi. 

Property (iii) is obvious, and property (iv) is a special case of (14. 9p . 
Finally, to check (v), we note that 

Fi{x;y; A) = ye{u^X,x/y;p)Rt^{X,ujx/y) = uxe{x/y,u'^X,u'^Xy/x;p). 
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By (14.31) . this gives 



[x; y) = uixu[x/y; p) - 



e{p; p 



2\ 



= i^rit; 1) - - ^dit,±pt;p') 



t 9{p\p^) 6{p;p'^ 

Finally, we use (14. 8 p to express Z^'^ in terms of 
Corollary 6.3. When ti are as in (14. 7p and n is even, then 

7^<^(f f f)- ,,n(n+l) ^(Au;^Aa;"+^p)^ 

^„ (t0,tl,t2j - -C^ ^(;^3.p3)n2+2„+3 



□ 



w/iz/e if n is odd 

^„ It0,tl,t2j - W ^(A3;p3)n2+2n+3 

,n^2. 2XK_ ^^,,2nxfl^ ^, ,n x 2. ^2m;^. 



X 6'(-Cc;"A'; p') lim ^. ".Y_^ - p— a;'"A^(-pa;"A'; p') lim 



7. Symmetric functions 

We are now faced with two problems: to construct a function satisfying all 
properties of Proposition 16.21 and to show that this function is unique. The 
key for solving both problems is to replace property (iii) by a generic vanishing 
condition. This leads to certain symmetric multivariable theta functions, which 
after uniformization become the symmetric polynomials (13.90 . 

7.1. Symmetric theta functions. We denote by Vn the space of BC^n+i theta 
functions of nome satisfying property (ii) of Proposition 16.21 For ti G C*/{t = 
p^t = t~^}, we denote by V^(ti, ■ ■ ■ ,tk) the subspace of Vn consisting of functions 
ijj such that, apart from the trivial zeroes at 1 and ±p, vanishes at ti, . . . ,tk 
(counted with multiplicity). Then, properties (i)-(iii) of Proposition 16.21 can be 
summarized as 

G K-i(WJ,P^_^. (7.1) 

n— 1 n— 1 

Lemma 7.1. The space Vn has dimension 2n + 1. The space Vn{ti, . . . ,tfc) has 
dimension at least 2n + 1 — k, with equality for k < n. 
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Proof. A basis for the BCsn+i theta functions of nome is given by 

ti-3n-l^^_p2i^6n+3. pl2n+6-^ _ ^3n+2-i^(^_p2i^-6n-3. pl2n+6)^ 1 < J < 3/2 + 1, 

see |RSct Proposition 6.1]. It is clear that the subspace V„ is spanned by the 2n + l 
basis vectors with j ^ mod 3. 

Since V^(ti, ■ ■ ■ ,tk) is obtained by imposing k hnear conditions on Vn, it has 
dimension at least 2n + l — k. For the final statement, note that by the quintuple 
product identity (14.41) . Vn contains all functions of the form 

When bi = ti for 1 < i < k and the remaining bi are generic, this function is in 
Vniti, . . . ,tk) \ Vniti, . . . , tfc+i). Thus, as long as we impose at most n vanishing 
conditions, each additional condition decreases the dimension by one. □ 

Although we do not need it in full generality, the reader may find the following 
characterization of the space Vn helpful. 

Lemma 7.2. Fix a G C*/{t = p^t = t~^} with a ^ —u. Then, Vn is the space of 
all functions ip that can be written 

^(t) = 9{t, ±pt; p^) [A 9i-t^, -cot^, at^;p'^)e{ait^^, ■ ■ ■ a„_it^^; /) 

+S^(M^^■■■M^';/)), 

with A,B E C and ai, . . . , ctn-i, 6i, . . . , 6„. G C*. 

Proof. It is easy to check that any BCsn+i theta function ip of nome p"^ can be 
written 

m = 9it, ±pt- p") {f{t') + tg{e) + t-^g{t-^)) , 
where / is a -D^+i theta function of nome p^, and g satisfies 

3(P% = ^^n9it)- (7.2) 

Since the term involving / has the desired form, we restrict to the case / = 0. By 
the quintuple product identity (14. 4p . G if and only if 

e{-pH-p')g{t-') =te{-pH-'-p^)g{t). 

Together with (17. 2p . this implies that g vanishes on all zeroes of ^(— t, —p'^t;p^). 
Factoring g(t) = t~^9{—t, —p'^t;p^)h(t) gives 

^(t) = t~^e{t, ±pt;py{-t^;p^) {e{-pH^;p^)+te{-pH-^;p^)) h{t^), 

where h is a Dn theta function of nome p^. 

Next, we observe that 6{—p'^t^]p^) + t9{—p'^t~'^\p^) vanishes on the zeroes of 
Thus 

0{-pH'; p') + e{-ph-'- /) = te{-t- p')kit), 
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where k is a D2 theta function of nome p^, that is, k{t) = C6{l3t^;p'^) for some C 
and (3 (depending on p). After simphfication, we conclude that 

^(t) = Ce{t,±pt;p^)e{-t^,-ut^,/3t^;p^)h{t^). 

Using (14. ip to write 

we arrive at an expression of the desired form. This shows that any can 
be expressed as indicated. The converse follows by similar arguments. □ 

In view of (17. ip . we are mainly interested in the space V^(ti, . . . , )f:2n)- Generi- 
cally, one expects it to be one-dimensional and spanned by the alternant 

^W = ,,.det (V',(t.)), 

with {ipjY^^ a basis of Vn and t2n+i = t- However, we have not found such 
expressions useful for our purposes. Instead, we will work with the following less 
symmetric determinants. 



^1(^1,^2,^3)= det {ijjiti)) 

l<i,j<3 



Theorem 7.3. Fixing a basis ipi, ip2, "^3 ofVi, define 

and, more generally, 
^'„(ti,...,t„,Mi,...,M„,z;) = — ^— ^ det ' 



e{y,±pv]p'^Y-^ i<Cj<n \u~'-^e{upf^;p^) 
Let (t/'^-"^)^"!^ be a basis ofVn- Then, 

^n{h,...,t2n+l)=C det J^Pfm, 

l<jj<2?i+l 

with C independent of each ti. In particular, is anti- symmetric in all 2n + 1 
variables. 

Note that is only defined up to a multiplicative constant, which we do not 
specify. 

Remark 7.4. If we choose the basis tjjj^^ as in the proof of Lemma [7. ![ then the 
alternant is a minor of the determinant corresponding to the BC^n+i Macdonald 
identity. This should mean that \E'„ can be interpreted as an affine Lie algebra 
character. Although this observation may have interesting consequences, we will 
not explore it here. In §7.4[ we will see that the trigonometric limit case p = 
corresponds to characters of the orthogonal and symplectic groups. The corre- 
sponding limit of Theorem 17.31 is closely related to some determinant identities of 
Okada P, see (Tml) . 

We divide the proof of Theorem 17.31 into a few lemmas. 
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Lemma 7.5. For any i = 1, . . . , 2n, the map ti i— )■ • • • , ^2n+i) is an element 

of Vn(ti, . . . ,ti, . . . , t2n+l)- 

The proof of Lemma 17.51 is straight-forward. It will follow from Theorem 17.31 
that the statement holds also for i = 2n + 1. 

Lemma 7.6. // the map t \l/„(t, ti, . . . , t2n) is not identically zero, then 

dim • • • ,^2n) = 1- 
Proof. Take to with \l/„(to, ti, . . . , t2n) 7^ 0. Consider the functions 

V'fc(t) = ^n(to, ti, • • • , tk-i,t, tk+1, . . . , t2n), k = 0,1, . . . ,2n - 1. 



By Lemma [73} these functions all belong to the space Vn{t2n)- Moreover, they 
are linearly independent since 

i'k(tj)j^O <^=^ j = k, j,k = 0, . . . ,2n - 1. 

Since, by Lemma [TTTl dimV"„(t2n) = 2n, we conclude that ipo, ■ ■ ■ ,4'2n-i span 

Vnit2n). 

Suppose now that / G Vn(ti, . . . ,t2n)- In particular, / G Vn{t2n), so we can 
expand 

2n-l 

fit) = J2 AkMt)- 

k=0 

The remaining vanishing conditions for / give Ai = ■ ■ ■ = A2n-i = 0. Thus, / is 
proportional to ipo, so dim . . . ,t2n) = 1- D 

Proof of Theorem I7.3[ Let 

. . . , t2n+l) 



C(ti, . . . , t2n+l) 



deti<ij<2n+l{lpj^\ti)) 



We first prove that C is independent of ti. By Lemma 17. 5[ the denominator 
and numerator are both in ^4(^2, • • • ,^2^+1)- If that space is one-dimensional, C 
is independent of ti. Otherwise, Lemma 17.61 gives C = 0. The same argument 
applies to the variables t2, . . . ,t2n- Thus, 

^„(ti,...,Wi) = C(Wi)^^.det {i^f\ti)), 

where it remains to show that C is independent of t2n+i- 
To complete the proof, we write 

C{v) '^riitl,...,tn,Ui,...,Un,v) 



C{ti) '^niv, t2,...,tn,Ui,..., M„, ti) 

and specialize m„ = wtn- It follows from the definition of that 



(7.3) 



\E'n(tl, . . . ,tn,Ui, . . . , Un-1, CUt„, v) 



UUtn'ufe{tltf,u^tf;p') 
0{v, ±pv; p"^) 



X -^^litn, Ujtn, t;)^„_i(ti, . . . , tn-l,Ui, Un-1, v) . 
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Using this in (17. 3p gives 



We claim that 



"^n-lih, tn-1, Ml, • • • , Un-1, v) 


^(ti, ±pti;p^)e{tltf; p^)^i{tn, OOtn, v) ^ 



e{v, ±pv- p2)^(t3^±3; cutn, h) 

To see this, consider the denominator and numerator as functions of v. They both 
belong to the space Vi(t„,a;t„), which is one-dimensional for generic t„. Thus, the 
quotient is independent of v and can be computed by letting v = ti. 
We have now reduced (17. 3p to 

C{v) _ . . . , tn-1, Ui,..., Un-1, v) 

C{ti) ^n-l{v, t2,..., tn-1, Ml, • • • , Un-1, h) ' 

By iteration, we conclude that 

C{v) ^ ^i(ti,m,M) ^ ^ 
C{ti) ^i{v,ui,ti) 

This completes the proof of Theorem 17.31 □ 

7.2. Uniformization. The next step is uniformization, see §4.11 We will work 
with the uniformizing map 

d{-pu;;pre{ujt^;p') 

' e{-uj;p^ye{piot^;p^) ^ ' 



and write 



Lemma 7.7. We have 



u^e{-l,-pu;p^) 
e{-p, -uj;p^) 



^(-1) = 1, ^(-a;) = C, 

C 



e(i) = 2c+i, a?) 



C + 2 

Proof. The first two identities are obvious. By (14. 6p . 

X ^ uj0ip,puj;p^)9{-pu;pye{st^;p'^) 

^^''^ ' se{-uj-p^YO{pus^,put^-p^) ■ ^ ' 

Plugging in s = 1 and t = —1, and using (14.21) . we obtain ^(1) — 1 = 2(^. The case 
s = p and t = — 1 similarly gives ^(p) — 1 = —'2C,{p)/C- Q 
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We denote by Wn the space of polynomials q of degree at most 3n such that 

m = e{t,±pt;py{piot^;py-qm), 

is an element of Vn- Then, the correspondence between ip and g is a bijection. We 
denote by Wn{xi, . . . , Xk) the subspace of Wn such that ip ^ Vniti, . . . ,tk), where 

Lemma 7.8. The space Wi is spanned by the three polynomials f , g and h defined 
m (133|) . 

Proof. By Lemma I7.2[ Vi is spanned by the functions 

e{ph^^;p^), e{-t^,-ut^,put^;p^), e{t^^-p^). 

By Lemma [7771 the uniformization of these functions are indeed, up to multiplica- 
tive constants, respectively 

(C + 2)x - C, {x- l){x - 0, x\x - (2C + 1)). 

□ 

Taking determinants of the basis elements fl3.6p . we get the following results. 
Corollary 7.9. For any elements qi, q2, qs ofWi, 

det {qj{xi)) = C{x2 - xi){x^ - Si)(x3 - X2)F{xi,X2,x^), 

l<i,j<3 

where C is a constant and F is as in (13 ■7p . 
Moreover, writing x = ^{t) and y = ^{u), 

u-^6{uH^^]p^)d{j)UJt^,puu^;p^) = C{y-x)G{x,y), 

where C is a constant and G is as in (13.81) . 

It follows that the function \E'„ from Theorem 17.31 is given by 

2n+l 

^n(ti,...,W) = C n 0{U,±pU;p')e{pu;tf;p^r U t-'0{tjtf;p') 

i=l l<i<j<2n+l 

X SniXi,...,X2n+l), (7.7) 

where C is a constant, Xi = ^(t,), and S'„ is the polynomial introduced in (13.91) . 
In particular, S'„ is symmetric in all variables. 
Finally, we note the following symmetry of Sn- 

Lemma 7.10. Indicating also the dependence on (, the polynomial Sn satisfies 

Sn{Xi , . . . , X2n+l] C ) = '72n^~n n Sn{xi, . . . , X2n+l', C) 

^1 ■ ■ ■ ^2n+l 
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Proof. This follows from an elementary computation, using 

Fix-\ y-\ z-'; C') = y, z- C), 

Q^x^y^ 

□ 

7.3. Recursions. Note that any minor of the determinant in fl3.9p is a determi- 
nant of the same type. Thus, any algebraic relation between minors (see e.g. [HI 
Chapter 6]) implies a relation involving the polynomials S'„. It is not our purpose 
to give an exhaustive list of such identities; we only mention a few examples that 
will be used below or otherwise seem of particular interest. 

First, we apply Jacobi's identity (sometimes called the Lewis Carroll formula 
after one of its proponents) 

x^n— 1,71 -y"n -yn—l \^7i 1^71—1 

where X = deti<ij<n{xij) and the other quantities are minors, upper and lower 
indices signifying omitted rows and columns. Choosing 

^ _ F{xi,yj,z) 
G{xi,yj) 

one obtains after relabelling the following recursion for the polynomials 5*^. 
Lemma 7.11. For x = {xi, . . . , X2n-i), 

(a -b){c- (i)5'„_i(x)S'„+i(a,6, c, c?,x) 

= G{a, d)G{b, c)S'„(a, c, x)S'„(6, d, x) — G{a, c)G{b, d)Sn{a, d, x)5'„(6, c, x). 

One of the Pliicker relations is 

IT" y — IT" y _L V" y — n 

^7i,n+l^n-l ^7i-l,n+l^n ^n-l,7i^n+l — 

valid for minors of an n x (n + 1) matrix. Specializing the matrix entries as above 
yields the following result. 

Lemma 7.12. For x = (xi, . . . , X2n-2), 

(c — d)G{a, b)Sn~i{b, x.)Sn{a, c, d, x) + (rf — b)G{a, c)5'„_i(c, x)S'„(a, b, d, x) 

+ {b — c)G{a, d)Sn-i{d, x)S'„(a, b, c, x) = 0. 

Combining Lemmas 17.111 and I7.12[ we obtain the following recursion. 
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Corollairy 7.13. For x = (xi, . . . , X2n-2), 

{b - d){y - a){y - c)S'„_i(a, x)S'„_i(c, x)5'„+i(?/, a, 6, c, d, x) 

+ (a - y)G{b, c)G{y, d)Sn-i{c, x)S'„(a, c, d, x)S'n(?/, a, 6, x) 
+ (y - c)G'(a, 6)5'„_i(a, x)S'„(a, b, c, x)S'„(?/, c, d, x) 

+ (c - a)G'(y, &)G'(l/, d)Sn{a, b, c, x)S'„(a, c, rf, x)S'„_i(y, x) = 0. 

Proof. The sum of the first two terms can be written 

{y - a)Sn-i{c, x){(6 - d){y - c)Sn-i{a,x.)Sn+i{y,a,b, c,d,x.) 

- G{b, c)G{y, d)Sn{a, c, d, x)5'„(?/, a, &, x) }. 
By Lemma [7.111 with (x, a, b, c, d) i— i- ((x, a), d, b, c, y), the factor in brackets equals 

-G(c, d)G{y, b)Sn{a, b, c, x)S'„(y, a, d, x). 
The sum of the last two terms is 

G{y, b)Sn{a, b, c, x) {(?/ - c)G{a, d)Sn-i{a, x)S'„(?/, c, rf, x) 

+ {c- a)G{y,d)Sn{a,c, d, x)S'„_i(?/, x)}. 
By Lemma 17.121 with (x, a, b, c, d) h> (x, d, a, y, c), the factor in brackets is 

{y - a)G{c, d)Sn-i{c, x)S'„(y, a, d, x). 
It is now clear that the sum of all four terms vanishes. □ 
Corollary 17.131 is particularly interesting in the case when 

X = (a, . . . , g, 6, . , 6 ), 

n— 1 n— 1 

c = b and d = a. As we will see, one special case is the three-term recursion for 
the polynomials P„. 

Corollary 7.14. For fixed a andb, the polynomials 

Sn{y) = Sn{y,a^_^^,b^_-2b) 

n n 

satisfy 

{b - a){y - a){y - b)Sn-i{a)Sn-i{b)Sn+i{y) 
+ ((a - y)G{b, b)G{y, a)5„_i(6)5„(a) 

+ {y- b)G{a, a)G{y, 6)5„_i(a)5„(6)) Sn{y) 
+ {b- a)Giy, a)Giy, 6)5„(a)5„(6)5„_i(y) = 0. 
We also consider a similar specialization of Lemma 17.111 
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Corollciry 7.15. For fixed a, b, and z, the polynomials 

Sn{x,y) = Sn{x,y,z, a,..^.,a ,b,..^.,b) 

n—l n—l 

satisfy 

(x - a){y - b)Sn-iia,b)Sn+i{x,y) 

= G{a, y)G{b, x)Sn{a, 6)S'„(x, y) - G{a, b)G{x, y)Sn{x, b)Sn{a, y). 

Two cases of Corollary 17.151 are of special interest. The first one is x = 6 and 
y = a, when it reads 

(a - bfSn-i{a, b)Sn+i{a, b) 

= G{a, bfSn{a, a)5'„(6, b) - G{a, a)G{b, b)Sn{a, bf . (7.8) 

This will yield Theorem 13.71 The second one is the limit case x a, y ^ b, where 
we first divide through with [x — a){y — b) and then use THopital's rule on the 
right-hand side. Let 

Tn{x, y) = Sn{x, ...,x,y,...,y,z). 



Since 



dT„, 



dx 



dSr, 



n- 



x=a,y=b 



dx 



x=a,y=b 



and similarly for y, the result can be expressed in terms of r„ and its derivatives. 
Corollary 7.16. The polynomials Tn = Tn{x,y) satisfy 

d^n on on 



(dG dG d'^G \ 1 / 
a M n H ( 
OX oy oxoy J \ 



dxdy dx dy 



Equivalent ly, if we let 



T 

J- n. 



then Tn satisfies the two-dimensional Toda molecule equation 

g2 



logr„ 



Tn-lT-n+l 



dxdy - 

We have recovered an instance of the well-known fact that this equation can be 
solved by determinants, which can be found already in Darboux's classic treatise 
[D| §378] (we owe this reference to Jacques Perk), see [HI ILSj for more recent 
accounts. 

The one-dimensional Toda equation 

, Tn-lTn+l 
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plays an important role in the analysis of the six-vertex model with domain wall 
boundary conditions [BFl IBLTI 16121 [BL3l iKZl [Zi]. This suggests that Corol- 
lary [7iT6] might be useful for proving Conjecture 13.141 

7.4. Trigonometric limit. In the trigonometric limit p — > 0, ( 17. 5p reduces to 
( = —2. We will show that the corresponding limit of the function Sn is related 
to symplectic and odd orthogonal characters. 

Let A = (Ai, . . . , A„) be a partition. Recall |FHt §24.2] that the characters of 
the Lie algebras 5p(2n) and so(2n + 1) are given by 

1 , /,-(A,+n-i+l) Xi+n-i+1 



,n , 1^ deti<i j<„ t • — t ■ 



-(Ai+n— i+i) ,Ai+n— j+i 



In this notation, we have the following result. 

Theorem 7.17. Let the variables ti, Ui, Xi and yi he related by 

Xi = -{l+ti + tT^), y^ = — TT- 

1 + Mj + 

Then, 

(yi Cy2m 



lim (C + 2)2'""-'"('"-i)^„ ( 

C^-2 \^ 



C + 2''"' C + 2 



2m 

2^{2m+n)n—m{m—l) 'j2mn—m(m—l) 



n 

+ + 

SD{4n-4m+3) /, , ^ 

^ ^{n-m.,n-m,,n-m.-l,n-m.-l,...,l,l,0)\''^^ • • • ' ''2n-2m+ly 

^ ^(m-l,L-l,. ..,1, 1,0,0) • • • 5 ""am), 



(7.9) 



lim (C + 2)(2"+i)"— f 

C^-2 \^ 



Xi, . . . , X2n— 2m; ^ _j_ 2 ' ' ' ' ' -\- 2 

2m+l 



2{2m+n+l)n-m22(2m+l)n-m2 TT ^ 

t\ {l + U^ + U-'r 
5o{4n— 4m+l) fj.j.\ 

(7.10) 

^ ^(n-m,n-m-l,n-m-l,...,l,l,0)v''l' ' ' ' ' ''2n-2m) 



X 



Bp(4m+2) / N, 

^(m,m-l,m-l,...,l,l,0,0)V^l' ' ' ' ' ^2m+lj- 



THREE-COLOUR MODEL WITH DOMAIN WALL BOUNDARY CONDITIONS 35 

Note that the correspondence between Xi and ti is the trigonometric hmit of the 
uniformization (I7.4p : 

\imi{ti) = -{l+U + tr^) =Xi. 

p->0 

For the variables Ui and Wj, we have instead 



C + 2 



1 + 2C 3 
^(mJ 1 + + 



which is a natural limit in view of Lemma 17.101 

We will use some determinant identities due to Okada [O]. Let 

W''{xi,...,Xn]ai,...,an) = det (x^"^ + ajX^^). 

l<i,j<n 

(In [O], ly" denotes the matrix rather than its determinant.) Then, 



det ' ^ _ 1 



i< 



ij<n\{l - Xiyj){yj - Xi)) n"j=i(l - ^iyj){yj - a^i) 

X 1/1, ...,?/„; ai, a„, 61, (7.11a) 

iW^{xuyj,z]ai,hj,c)\ _ (1 + c)"-i 



i<i,i<n V (1 - a;i%)(2/i - ^i) J n"i=i(i - Xiyj){y. 



X ly ''"^ (xi, x„, yi, ...,?/„, ai, ... , a„, 61, ... , 6„, c). (7.11b) 

Okada used these identities to enumerate certain symmetry classes of alternating 
sign matrices. 

The characters appearing in Theorem 17.171 can be obtained as special cases of 
Okada's determinants. 

Lemma 7.18. The following identities hold: 

iy^"+^(ti, . . . , t|„^i; — ti, . . . , — t2n+l) _ B0(4n+3) /, , n 

n?:r tr(i - u) i{,._.,,^..ut, - t.)(i - t.t,) " ^("— -.■■■.^.i.o)(*i^ • • • ^ wj, 



:-i)"iy^"(tf,...,tL;-ti,...,-t 

n?:iir'(i - i.') ni<.<,<2n(^. - m - uh) 



2n ,n-l(. ,2Nrr /'+ +\^^ ++\ '^(n-l,n-l,...,l,l,0,0)'^''l'---'''2n;, 



Vr^""'"^(ti, . . . , t\^j^i] tf, . . . , t^n+i) _ sp(4n+2) /, , n 

n£rt«(i - 1?) ni,.<.,2..i(^. - m - ut,) - ^(-"— .-.maoI*!' • • • ^ wi. 

( — l)"iy^"(tf , . . . , —if, . . . , —tin) _ so{4n+l) / , , x 

A(n,n-l,n-l,. . .,1,1,0) v''!' " " ' ' ^2nJ- 



n?:i - ni<.<K2n(^. - - t.^.o 



<i<'j<2n^ 
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Proof. This is straight-forward, and we only provide some details for the first 
identity. By definition, 



2n+l 

- - 6n + l / \ 

-L-L ' l<i,j<2n+l V * ' / 



=1 

where //j — 3n — 3j + 7 /2. In order to make all /ij positive, we multiply the last 
n rows by —1. To sort the Hj in descending order, each of those rows must be 
moved an odd number of positions upwards. Taken together, these two operations 
do not change the determinant, and we arrive at the desired result. □ 

We will need the following elementary identities. 

Lemma 7.19. In the notation above, the following identities hold: 

{l-tlt2){t2-ti) 



X2 - Xi 



y-2-yi 



tlt2 

3(1 - Ui){l - U2){1 - UIU2){U2 - Ui) 

{l-u\){l-ul) 
2{tl-t\){l-t\tl) 



lim iC + 2?r fr \ - 



2 



72(1 - UiYjl - Uifjul - ul){l - ulul) 
C'A"V^ '^'^VC + 2'C + 2y~ {l-ulf{l-uiyiu2-ui){l-uiu2) ' 

2W%tltltl;-ti,-t2,-ts) 



c->-2^ ^ vc+2'c+2; 



lim F{xi,X2,X3) 



' ' nlMi-ti)Ui<i<j<sitj-ti)ii-titjy 

lim (C + 2)F (x X ^\ = l2u,{l-u,)W\tltl-tl-tl) 
^nm^^C + yx^,X2. ^^^j ^3^^^^^ _ ^^^^^^^ _ ^^^^^^ _ ^^^^^ _ ^^^^^ , 

Cvi Cy2_\ 

72 uiU2W^{ul, ul; -Ui, -U2) 



(1 + Ui){l - Ul){l + U2){1 - UI){U2 - Ml)(l - U1U2) ' 

2 r / Cvi Cy2 Cys 



lim (C + 2YF . ^ 

C->-2^^ ' VC + 2'C + 2'C + 2, 

72W^3(^3^^3^^3._^2^_^2^_^2) 
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Proof of Theorem 17.171 To prove fl7.9l) . we relabel the variables, considering in- 
stead 



lim (C + 2)2™"-™('"-i) 



x5„ 



Applying fl3.9p leads to the block determinant 



Xr 



• (7.12) 



det 



( -P ( ^ <yj_ y 

\^C+2' C+2' 
F I T- ^'^^ 7 



v 



Gix- ^ 



J 



By Lemma [7.191 the off-diagonal block;s vanish in the limit C — )■ —2, so the deter- 
minant splits as the product of the diagonal blocks. It follows that f l7.12p can be 
written 



lim TT 

-^-9 J-J- 



C^-2 



m+l<j 



x(C + 2r(-+^)5„(^, 



C + 2'C + 2' 

,yn,z). 



C + 2' 



Computing the prefactor using Lemma 17.191 we are reduced to proving the two 
cases m = and m = n of (17.91) . that is, 

^hm^S'„(Xi, . . .,X2n+l) = 2'' X(^^n?7t^l,n-l,..,l,l,0)(^l' " " • '^2n+l), 



^l_im^(C + 2) ^"^^ + 2'---'C + 2''^ 



2n ^ 

2n(2n+l)on.(n+l) TT t «P(4n) / ^ 



This is now straight-forward, using Lemma I7.19[ Okada's identities (17. lip and 
finally Lemma [7.181 
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The proof of (17. 101) is similar. We consider the limit 



hm (^ + 2)(2"+i)"-"^' 



X Sr: 



',y Vm+l) • • • ; Vny 



Incorporating a factor {( + 2)™ into the determinant, we consider 

/(C + 2)F(|f|,g|,^) (C + 2)F(g|,y,,^U 



hm det 



v 



^ \C+2-> C+2 



G'(a;i,|/j) 



In the hmit, the upper right block vanishes. Similarly as before, this reduces the 
proof of f lT.lOp to a straight-forward verification of the special cases 



lim (C + 2)"5'„ ( xi, . . . , X2n, 



nn{n+l)Qn 

^ "J SD(4n 



'''-(n,n-l,n-l,...,l,l,0)^ 



'1 + Ml + 



L) • • • 5 ^2n)i 



^ito(c+2r''.«)s„(^,.... 



2?i+l 



C + 2 



2n+l . 

9n(2n+l)on(n+l) TT ^ sp(4n+2) / x 

11 (1 + Mi + U-^Y ^(n,n-l,„-l,...,l,l,0,0)l«l' • • • ' «2n+lj. 



□ 



As a consequence of Theorem I7.17[ we have the following identities. 
Corollary 7.20. One has 



lim 1 + - 

C->-2 V 2 



. „2 r (n-l)^ 1 

CA ^L^J o , . , . ^ ^ 

C + 2'---'C + 2 



^J2C + 1,...,2C + 1, 
V ' 



ri+l 



2" 3" An+i, n even, 
2'^'3'^a+i, n odd. 
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lim 1 + ^ ^J2C+1,...,2C + 1,-^,...,-^ 
C->-2 V 2/ v^^ ^ + 2 C + 2 

n . ' 

71+1 

2_-i n(n-2) 

2" 4 Cn+i, n even 

_2"'-i3^^A„+i, n odd, 

where An and Cn are as in (II. ip and ( 12.6P . 

To verify these identities, we note that the hmits correspond to special cases of 
Theorem 1 7 . 1 71 where ti = Ui = 1 for all i. Recall that, if fii = Xi + n — i, then |FH] 

sp(2n)/, ,N _ Yll<i<j<nil^i ~ /^i) Yll<i<j<nil^i + f^j + 2) 



xr^^(i,...,i) 



2"l!3!---(2n-l)! 



so{2n+l)/i ^ N _ Y[l<i<j<ni(^i l^j) Y[l<i<j<nil^i ~^ f^j ~^ ^) 

I '■••'J- l!3!---(2n-l)! ' 

The relevant special cases can be simplified as 

SD(4n+3) 1 ^ = TT (^^+ 

A(n,n,n-l,n-l,.,l,l,0)'.-^' • • • ' -^^^ 1 1 (2n + 2k + iV 

sp{4n) , X _ on(n-l) TT (6^ + 4)! 

X(n-i,n-i,...,i,i,o,o)l^'---'^J 11 (2n + 2fc + 2)!' 

fc=0 ^ ' 

A(n,n-l,n-l,..., 1,1,0,0)1-*-' " " " ' -"-J 22"+! H {2n + 2k + 1)! ' 
^so{4n+l) .X <^2non(n-l) TT (6^ + 5) (6/i: + 3) ! 

X(n,n-i,n-i,...,i,i,o)l^'---'^^ ^ 11 (2n + 2A; + 2)! ' 

Using these identities, it is straight-forward to derive Corollary 17.201 from Theo- 
rem 17.171 

8. The polynomials P„ and p„ 

In this Section, we show how specialization of the variables in 5*^ lead to the 
polynomials P„ and and deduce a number of properties of the latter two 
systems. 

8.1. Elementary factors. Our first task is to identify some elementary factors 
that appear when specializing the polynomials Sn- We will use the notation 

,2" 



5n 



n 



4 

cf. (13. 2p . For later reference, we mention the identities 

5n+i + 5n-i - 2(5„ = odd), (8.1) 
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^^^^^1 =n'- 25„-i = 25„ + n. (8.2) 
Proposition 8.1. There exist polynomials Pn{x,(), PniC) o,'>T'd UniC) such that 

c c 



(^1 + ^^(x, 2C + l,...,2C + l , 



^2' 'C + 2 



n 



(l + gn( 2C + l,...,2C+l , 

n+1 



C + 2'---'C + 2 



n 



(-i)("^')c'(i + cr'(i + 2C)^"(i + ^] 'p„(c), (8.3b) 



Sn- 



and, for n>l, 



c 



n(n— 1) 



1 + >gn( 2C + l,...,2C + L 

n+2 



'C + 2''"'C + 2 



n-l 



(_l)("^^)2mC"'-i(l + C)"'(l + 20'"+^ (1 + M " i/„(C). (8.3c) 



Applying Lemma [7. 101 to fl8.3ap gives, using also ( 18. 2p . 

P„(a;,C) = x"c['^]Pn(l/x,l/C) 
Similarly, applying Lemma [7. 101 to (]8.3bp gives 

/ >\n(ri+l) 

U + y ^nj 2C + l,..^.,2C + L 



•4) 



"'C + 2''"' C + 2 



n+1 



(-i)m2-[^]c"("+i)(i + c)"'(i + 20^"- f 1 + ^ ] ^"(0, 



where 



We will see in Proposition 18.91 that degp„ = n(n + l)/2, so the notation agrees 
with ( 13. 3p . Moreover, we clearly have 

P42C + 1,C) = (1 + 2C)[5W(C), (8.5a) 
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Pn c) = c(c + 2)-mp„(c). (8.5b) 

To prove Proposition 18. ![ we first observe that S'„(xi, . . . , X2n+i) is a polynomial 
of degree at most n in each Xi. This is obvious for the variable z in f l3.9p . and 
thus holds in general by symmetry. It follows that the left-hand sides in (18. 3 p are 
polynomials in x and C,. We need to show that they vanish of appropriate degree 
at the points C = 0, —1, —1/2 and —2. 

We first consider the function 

^^(l + Ca;i, . . . , 1 + Ca^n, Cz/i, • • • , Cz/n, z). (8.6) 

Expressing it as in fl3.9p . we observe that since | G{1 + (x,(y), the prefactor 
nrj=i G{1 + (xi, Hj) is divisible by C^"" . The product 

n {{l + Cx^)-{l + Cxi)){Cyj-Cy,) 

l<i<j<n 

contributes a factor ^"■("■^i) to the denominator. Finally, since ( \ F{1 + (x, (y, z), 
each matrix entry has a single pole at C = 0, so the determinant has a pole of degree 
at most n. In total, we conclude that f l8.6p is divisible by -"■("■-!)-" = (^^ . This 
implies the same statement for the left-hand sides of f l8.3ap and (I8.3bp . 
In the case of fl8.3cp . we consider instead 

Sn (^1 + C^^l, . . . , 1 + CXn-2, Xn-1, Xn, CVl, C?/n-2, Vn-l, " (^•'^) 

Using that ( \ G{x, (y) and that 

C \ 2C2(C + 1)^ 



G X 



C + 2J (C + 2)2 

for all X, obvious modifications of the previous argument shows that fl8.7p . and 
hence also the left-hand side of fl8.3cp . is divisible by 
The case C = — 1 is simple. Considering 

^„(-l + (C + l)xi, . . . , -1 + (C+ l)xn, -1 + (C+ 1)1/1, . . . , -1 + (C+ l)Z/n, z), (8.9) 

and observing that 

(C + 1) I F( - 1 + (C + l)x, -! + {( + l)y, z) , 

(C + l)'|G(-l + (C + l)x,-l + (C + l)z/), 

one checks that (18. 9 p is divisible by {( + 1)" . This implies the corresponding 
statement for the three left-hand sides of (18.30 . 

The vanishing conditions at ( = —1/2 and ( = —2 follow from the following 
Lemma. 
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Lemma 8.2. As a polynomial in (, 

Sn{{2( + l)xi, . . . , (2C + l)Xk,Xk+l, . . .,X2n+l) 

is divisible by {2( + lY'^-'^, and 
zs divisible by (C + 2)'^'=-i. 

Proof. The second statement follows from Theorem 17.171 (note that we only need 
the very first step in the proof, not the formulas involving Lie algebra characters). 
The first statement then follows using Lemma 17.101 □ 

8.2. Recursions. We will now specialize the recursions for Sn given in §7.31 to 
obtain recursions for P„ and Pn- 

If we let a = 2C + 1 and 6 = C/(C + 2) in ([7S]), using (^Ml and the identities 

G(2C + l,x) = 2(C + l)V, 

2C + l-^ = ^^^ 
we obtain the following result after simplification. 
Proposition 8.3. The polynomials pn and yn satisfy 

(C + i)V+i(CK-i(C) = C^'PniOyniO 

/ A\ 1+X(" odd) 

\1+X("cvcn) (ill) „ (/■\2 



+ (1 + 20^+^^^"^™"^ (^1 + IJ PniCr- (8.10) 

Equivalently, (13. 5p holds, where Y = yn- 

By induction, we obtain the following consequence, which is a key result in our 
approach. 

Corollary 8.4. The polynomial pn satisfies Pn(0) = 1; in particular, it does not 
vanish identically. 

From this, another key result follows. 

Corollary 8.5. The space Vn{l, . . . , l, p, . ,p ) defined in §7.11 has dimension 1. 

n n 

Proof. By Lemma [7.61 and (17. 7p . it is enough to show that the polynomial 

c c 



5„(a;,2C + l,...,2C + l, 



C + 2''"'C + 2 



does not vanish identically. By ( ]8.3bp , this follows from Corollary 18.41 □ 
We remark that Corollary 18. 51 can alternatively be deduced from Theorem 17.171 
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Proposition 8.6. The polynomials Pn{x,Q O'f^ determined from the starting val- 
ues Pq = 1 and Pi = X + C, by the recursion 



Proof. Letting a = 2C + 1 and 6 = C/(C + 2) in Corollary [TH gives flSTTTl) after 
simplification. The polynomials are determined by this recursion, provided that 
An{x, C) does not vanish identically. By (18. 5p . this follows from Corollary 18.41 □ 

At this point, a remark on our approach is in order. It is easy to see that 
if a system of polynomials satisfies (18. lip , then it can be used to construct (via 
uniformization) functions satisfying the properties of Proposition 16.21 The 
reader may ask why we do not take (18. lip as our starting point, thus avoiding 
introducing the multivariable polynomials Sn- The problem is that, to use (18. lip as 
a definition, we must know that An{x, () does not vanish identically. Moreover, we 
need to know that is uniquely determined by the properties of Proposition 16.21 
Both these facts have been deduced from Corollary 18.41 If one could find a way to 
deduce inductively from (18. lip that P„(2C + 1,C) and Pn(C/(C + 2), C) never vanish 
identically, then the proofs of our main results could be significantly shortened. 

8.3. Further properties. Using Proposition 18. 6^ we can obtain quite detailed 
information on the polynomials P„. 

Proposition 8.7. The polynomial Pn{x , () co,n be written 



A„(X, C)Pn+l(x, C) = B^{X, C)P„(X, C) + Cn{x, C)Pn-l(x, C), 



(8.11) 



where 




Cn{x,0 




n 




(8.12) 



k=0 



where is a polynomial of degree 2Sn, with leading coefficient 




(8.13) 
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In the notation f l3.3p . these polynomials satisfy 

fkiC) = /:-.(C). (8.14) 
For k> [{n + l)/2], fJ!{C) is divisible by (C + 2)^~[^] . Moreover, 

[~] 

/:(C)= (1 + ' Pn-i(C), (8.15a) 

/:-i(C) = 2 + Ij iC + n + l) p„„i(C). (8.15b) 

Further results follow by applying (I8.14p to the other statements. In particular, 
from fl8.15ap we obtain 

PniOX) = C f^^V'^P„_i(C). (8.16) 



Proof. From the proof of Proposition 18. H and also from Proposition 18. 6[ it is clear 
that P„(a;, Q is of degree at most n in x, so we can write it as in f l8.12p . where a 
priori /^(C)C"~'^ is a polynomial. 

To show that is a polynomial, we consider the function 

F^{x) = lim = y lim fj:{C)x\ 

^ fc=0 

We will prove by induction that the limit exists and equals 

in standard hypergeometric notation. To this end, we divide (18. lip by (^^" and 
then let C 0. We need that, by f l8.5ap and Corollary 18. 4 [ 

P„(l,0)=p„(0) = l, 

and we also write 

,_^^P,.(C/(C + 2).C)^ 

Then, the resulting identity simplifies to 

(a;-0F„_i(l/2)F„+i(x) 

= (^x'F„_i(l/2) +{^-^ ^n(l/2)^ F„(x) - x'F„(l/2)F„_ii 
We must show that this recursion is solved by (18.170 . Equivalently, since 



x). 



2" V J 2"+i V ^ + 1 
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we need to check that 



1\ , , ( ^ 2n + l f l\\ , , 2n + l ^ 



f:-M = -j : ). (8.1^ 



X--] Fn+lix) = [x' + ~ 2; J ^"^"^^ ~ T^TT ""^^"^^ 

which can be done by a straight-forward computation. 
The equation fl8.17p shows that is a polynomial, with 

1 /n + P 

It follows from (18 ■4p . using also (I8.2p . that 

mo = e'-Jn-kivo- 

Together with f l8.18p . this proves that fj! has degree 2Sn, as well as fl8.13p and 
The statement on divisibility by powers of C + 2 is equivalent to saying that 



c- 

exists finitely. This is a consequence of (18 .Sap and Theorem 17. 171 

To prove fl8.15p . we pick out the coefficients of x""*"^ and x""*"^ on both sides of 
(18.1ip . On the right, neither C„ nor the second of the two terms in Bn contribute. 
Using (18.51) and simplifying, we find that 



Pn-i(C) x'^ - 



c 



C + 2 

X{n, odd) 



+ 2c + x^ (/:i'(c)x"+^ + c/r'(c)a:") 



1 + ^ ) Pn(c)x'(x - 2c - 1) imox'' + c/„"-i(c)a:"-^) + ■ 



2, 

where the ellipsis denotes terms of order at most n + 1 in a;. Picking out the top 
coefficient gives the recursion 

/ /-^Xin odd) 

which is solved by (I8.15ap . Picking out the next coefficient then yields 

/ X(n odd) / [^] 

Pn-liOf:^\0 = (1 + Ij Pn(C)/:-l(C) + 2 + ij Pn-l(CK(C), 

which is solved by fl8.15bp . □ 

For the formulation of the next result, it is convenient to introduce the poly- 
nomials 

6(r)- F f'n/2,-in-l)/2 
Mx)-2F,\^ n + 3/2 '"^ 

/ . . p f-n/2, -(n + l)/2, (n + 5)/4 
Mx)=3F,l n + 3/2, {n + l)/A 
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of degree [(^ + l)/2], respectively. 

We will need the identities 

(3a; + l)Mx) = l^^'^.V/o ""^.t^ + (x - ifK-iix), (8.19) 
[2n + l)[2n + 6) 

(3X + irUx) = + - l)'^n-l(x), (8.20) 

(2n + l)(2n + 3) 

^{2n + l){2n + 3) ^'^"+1^^) = ~ ^^'^^^^^ + " 1) Vn-i(x), (8.21) 
which are straight-forward to verify. 

Proposition 8.8. For special values of(, Pn{x,() nT-dy be expressed as 

P„(x,0) = x", (8.22) 

Pn{x,-1) = (_l)x(n=2mod4)25„-i(3,_ ^g 23) 



P„(x, 1) = 2^— 3[t] A„+,(l + a:)>„ [^^^^^^^^ 



Pn{x, -2) = <^ 



/ a; + 3 \ 

(-3/4)2y4„+i(l - x)2 0„ I -y 1 , n et;en, 

-2-'^-i3"^a+i(l-x)^^„f-^^y norfd 



(8.24) 



(8.25) 



,3(1 -x) 

Applying fl8.4p to fl8.25p one may also evaluate P„(x, —1/2). 

Proof. The identity ( KT2^ follows from f lSTT^ and f lSTTB]) . 
If we let = — 1 in fl8.1ip . we obtain after simplification 

(_^)x(n even)2X(n °dd)p^_^(_^^ -l)P„+l(x, -1) = (x - 1)P„(-1, -l)P„(x, -1). 

It is easy to check that this is solved by fl8.23p . using (18.1 p to simplify the exponent 
of 2. 

Similarly, if we let C = —1 in (18. lip and substitute (I8.24p . we are reduced to 
the recursion (I8.2ip . For this computation, it is useful to note that 

AnAn+2 _ 3(3n + 2)(3n + 4) 
^2^1 ~ 4(2n+ l)(2n + 3) 

and that 

_(3x-l)(x-3) _ 16x _8(x2-4x + l) 

^~ 3(l + x)2 ^ y-'^- -3(3; + 1)2' (a; +1)2 • 

To prove IK25\\ . we muhiply f lHUj) by (1 + C/2)['^] and then let C ^ -2. This 
gives 

^_3)x(n cvcn)(^ ^ 3)X„_iP„_i (-3, -2)P„+i(x, -2) 

= X„ (-9P„_i(-3, -2)P„(x, -2) + x2p„(-3, -2)P„_i(x, -2)) , (8.26) 
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where 

It follows from (18 .Sap and Corollary 17.201 that 

l)52-a^„ neven, 
^ -2-M„+i, nodd. ^ ^ 

It is now straight-forward to check that substituting (18.251) in (I8.26P yields (I8.19P 
when n is odd and (I8.20p when n is even. For this computation, one needs the 
identities 

An+2Cn ^ 3(3n + l)(3n + 4) 
A„+iC„+i 4(2n+l)(2n + 3)' 
A„C„+2 _ 3(3n + 2)(3n + 5) 
~ 4(2n+l)(2n + 3)' 

and also that, 

x + 3 4x 4 
2^ — 1 = 777:; r, 3z + 1 



3(l-a;) 3(l-a;)' 

□ 

Proposition 8.9. The 'polynomial "Pn has degree n{n + l)/2 and leading coefficient 
2 J ( j . Moreover, pn assumes the special values 

J9„(-1) = (_l)x(n=l mod 4)2<5n+i ^8 28) 

p„(l) = 2^"+M„+i (8.29) 

/ )^n+u n even, 

l^(-l) 2 6„+i, n odd, 

l(— 1) 2 2 '2^'^+'-) An+i, nodd. 

Proof. The first statement follows (I8.13P and (jHTTSaj). The evaluations (I8.28p . 
flOg]) and dOU]) follow from Proposition ESI using (K5^ . Finally, by (l83b|) . 
fICT]) is equivalent to (KTHf . □ 

Writing p„(C) = XlfcST^''^^ we have simple formulas for oq (Corollary 18. 4 p 
and an(n+i)/2- It is easy to deduce from Proposition 18.31 that 

7n{n + 2) 



n even, 



ai 



7(n2 + 2n + 3) 

, n odd. 
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It also seems that 

n\7n+10) f2n + 2 



Q n(n+1) -I 
2 ^ 



2("+8)/2(n + 2) + 1 

(n + l)(7n2 + 3n-6) /2n + 2 
2('^+7)/2(n + 2) Vn + 1 



n even, 



n odd. 



One can probably prove this using Proposition I8.6[ though we have not worked 
out the details. It does not seem that the other coefficients admit such simple 
expressions. 

Finally, we use the ffist part of Proposition 18.91 to pick out the leading term on 
both sides of f lS.lUp . This allows us to compute the leading term in the polynomials 
Un- In particular, we may conclude that degY = degA in (13. 5p . 

Corollary 8.10. The polynomial i/n has degree "'-"''''^ + 2 and leading coefficient 

{2n + 2)\(2n)\ 



2n+x(n odd)(^^ + 2)!(n + l)!2n! ' 

8.4. Zeroes. In this Section we prove Propositions 13.12] and [STTSl 

Proof of Proposition I3.12[ We proceed by induction on n, treating the cases of 
odd and even n separately. Assume that the statement holds for the zeroes of p^n 
and P2n+i- Let ai, . . . , a„ denote the real zeroes of p2n^ bi, . . . , bn+i the real zeroes 
of P2n+i, and Ci, . . . , c„+i the real zeroes of P2n+2- We assume that 

6i < fli < 62 < 02 < ■ ■ ■ < a„ < bn+i < -2, (8.32) 

and need to prove 

61 < Ci < 62 < C2 < ■ ■ ■ < bn+l < Cn+1 < "2. 

Since, by Proposition 18. 9 [ p2n has positive leading coefficient, it follows from (I8.32p 
that {—l)^~^^~^^P2nibi) > 0. Moreover, by Proposition 18. 3[ if P2n+i(C) = 0, then 

(C + 1) W(C)P2n+2(C) = (1 + 20 (l + P2n+l(0'. 

We conclude that (— l)"^V2n+2(^i) > 0. Thus, p2n+2 has one zero between each 
consecutive pair of points 6j, Since, by (I8.30p . P2n+2(~2) > 0, the remaining 
zero lies between bn+i and —2. The induction step for odd n is similar, and we do 
not give the details. □ 

Lemma 8.11. Fix ( in the interval —2<(<—l/2 and assume Conjecture 13. Ill 
Then, Pn{() is negative if n = 1 mod 4 and positive else. Further, Pn{C) is 
negative if n = 2 mod 4 and positive else. Moreover, P„(0,C) is positive if n = 
mod 4 and negative else. Finally, the leading coefficient of Pn{x,() is negative 
if n = 2 mod 4 and positive else. 
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Proof. Since we assume that pn does not vanish in the interval —2 < ( < —1/2, 
PniO has the same sign as 1), so the first statement is obtained from f l8.28p . 
The sign of pn{C) is determined as a consequence. The remaining statements follow 
using f l8J[5al) and f l816|) . □ 

Proof of Proposition 13.131 For fixed (, with —2 < C < — 1/2 and ( ^ —1, let 
ai, . . . , Qn-i denote the zeroes of Pn~i, bi, . . . ,bn the zeroes of P„, and ci, . . . , Cn+i 
the zeroes of Pn+i- By induction, we assume 

< 6i < ai < 62 < 02 < ■ ■ ■ < ttn-i < bn, (8.33) 

and prove 

< Ci < 61 < C2 < ■ ■ ■ < 6„ < Cn+l. 

Using Lemma IS. lit we deduce from f l8.33p that 

(_iy+X(n-l mod 4)p^_^(^^^^)>0. 

It follows from (IHTTj) that, if P„(x, C) = 0, then 

"^'^ ' " (X - 2C - i)((C + 2)x - C)p„-i(CK-i(C) ""'^ ' 

Again using Lemma [8. IH 



> 0. 



Pn-l(CK-l(C) 

Combining these facts we find that 

Thus, Pn+i has a zero between any two consecutive bi. Moreover, yet again using 
Lemma [8.111 P„+i(0, C)-Pn+i(&i, C) < 0, so there is an additional zero between 
and 61. Finally, P„+i(6„, () and the leading coefficient of Pn+i{x, () have opposite 
signs, so the final zero is to the right of 6„. □ 



9. Return to three-colour model 

9.1. Uniformization of $„. After the long detour in ^ and ^ we are now ready 
to apply our results to the three-colour model. First, we express the function 
introduced in ^ in terms of the polynomials P„_i. The following identities will 
be useful. 
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Lemma 9.1. With ( as in (17.51) . 



2C + 1 
C + 2 = p 

c + i = - 



9i-l,-oj;p')9icu;p^y 
0{—p, —poj; p'^)6{j)uj] p'^y 

e{p, -puj]p^) 



C-1 



6{-p,puj]p'^y 

6{p,puj]p^)d{u]p'^Y 



(9.1) 
(9.2) 



J(— p, —pu)] p'^)9{—u; p^Y 

Proof. All four identities follow from Lemma 17.71 in the last two cases by writing 
C + 1 = - e(-c^), C - 1 = a-^) - and using (HD. □ 

We note in passing that Lemma [9.11 can be used to prove (15. 5p . It is enough to 
check that, after elementary simplification, 



(2C + 1)(C + 2)(C-1) 



2\4 



pu)- 



6{u;p^y^6{pu;p 



{p; p) 



this implies (15. 5p in view of (I3.13p . 
Lemma 9.2. One has 

c(c+lr = 2r^^ 

where 

9{p] p^)6{—puj;p'^y 
^ 6{pu;p'^)6{—p;p'^y 

Proof. This follows easily from (19. ip . using (14. 2p . 

Proposition 9.3. In the notation above, 

$„(t) = Bj{t,±pt-py{t^,pt^,put^;p'r-'p^_Mt)x), 

where 

f 1 

pu 



□ 



(9.3) 



5. 



—p\p^)Q{j)UJ\p^Y"^ "^1] 4 (2^ + 1) 2 



2 ' 



[ e{p-p'^)e{pu-p^Y^'^7]'^{2C + 1)"^ ' 



n even, 



n odd. 



Proof. Both sides of (19. 3 p belong to the space 

n— 1 n— 1 



which is one- dimensional by Corollary 18.51 Thus, (19. 3p holds for some constant 
which we compute using part (iv) of Proposition 16.21 We have 

= i-ir+'B.,co'--9{-pu; p')d{p- p'r-'9{u;- p'f^-'dipuj- p'f^-'P,,^,{0, (), 
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$n-l(t) 



2n-3 



t=p 



By (l83b|) and (l8A6|) . 



^'n-l(0,C) 



P„_2(C/(C + 2),C) 



(C + 2)(2C + 1) 



We conclude that 

Br, 



Bn-l 



2\2n-4 



e{-pu;p^)e{p,pu;p^r-^ C V(C + 2)(2C + 1) 



6{—p; p'^)6{pu; p'^Yr] 2 



n even, 



n odd, 



I e{-uj]p^)e{puj]p^yr]—{2c + i) 

where we used Lemma [OTT] to simphfy the expression. It is clear that this recursion 
can be solved as indicated; the starting value Bi = ijj/9{p;p'^) follows from part 
(v) of Proposition I6.2[ □ 

Corollary 9.4. When ti are as in (14.71) . 

z„ (to,^i,^2)- ^2 ; 

^(p;p2)ry-0(A3;p3)r^2+2n+3 

X (ctp„_i(C)e(-pa;,-a;"A2;p2) _^i-Ap„_i(C)e(-a;,-pa;"A2;p2)) 
for n even, while for odd n 

ryscr N_ ^(Aw^Aw"+i;p)2 

Z„ (to, ^1,^2)- ; 

e{p; p2)^ — ^(A3; p3)n2+2n+3 

X (p„_i(C)^(-p,-u;"A2;j92) -u;-"AC^Pn-i(C)^(-l,-pc^"A2;/)) • 
Proof. In Corollary 16. 3 [ express as in Theorem 19.31 Using (18. 5p . we see that 

where 



D„ = {-ir+'e{-p;py{pio;p'y^-'{2C + l)[—\B^, 
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To complete the proof, we use Lemma 19.11 to write these expressions as 



pu) 2" 5 even 



Dr 



9{p;p'^)r]- 
e{-p;p'^' 



e{p;p'^)r]' 



— , n odd, 



9{-puj;p'^)C"^ 



n even 



pg(l-n)-^.-n(n+l) ^(-l;P^)C ^ 



2—;-, odd. 



4 



□ 



Using fl8.30p and (18.311) . one may check that in the trigonometric case p = 
(which imphes r] = —1 and C = —2), Corollary 19.41 reduces to (12. 5p . 
When n = 2, we obtain the following important consequence. 

Corollary 9.5. In the notation of Lemma [5. Ij, 

C2 + 4C + 1 



t{p) 



V 



By Lemma [9. 2 [ it follows that, when tj are given by (14.71) . 



T(t t ,. ^ic'+^c+iy 



(9.4) 



Proof. Noting that 



d{\;pf + d{\uj^;pf 
e{X;pye{Xuj;py^e{Xu;^;py 



and that 



, , ^^(A; pf + 9i\Lo; pf + e{\u^-pf ^(A; pf + e{\uj^-p) 

T[P) = TTTT^ = lim 



^(A3;p3 



^(A3;p3 



we express r(p) using the case n = 2 of Corollary 19.41 After simplification, we 
obtain 



6{uj]pYd{p]p^)r] ^'^ ^'^ ^ ^ ( C + )) 

C'(C + 3) - (3C + 1) = (C - i)(C' + 4C + 1), 

applying (19. 2p and using (I4.2bp . we arrive at the stated result. 



Factoring 



□ 
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9.2. Return to the polynomials g„ and r^. Combining Theorem 13 . 1 1 and Corol- 
lary l9.4t one easily recovers two of our main results: Theorem 13.21 and Theorem l3.6[ 
As a starting point, we note the theta function identities 

— ^7 , (9.5a) 



X9{-pu''X';p') 

9i-pu; p'f9iXu'^^; pf - 9{-p- p^)9{\uj^-+\ Xuj'^^+'; p) 

9i^' ■ ^ ^ 

These can be obtained as special cases of (14. ip . or simply by noting that since 
they relate three functions in a two-dimensional space, it is enough to verify them 
at the points A = cu" and A = w'^"'"^. 

Using (19. 5p in Corollary 19.41 we obtain an expression of the form (15. ip . where 
for n even. 



^ bn-l(C) - C^Pn-l(C)) 

9{p;p'^)9{u;pyri * 

Pn-l(C) -C^Pn-l(C) 



(1-C 

u^9{-u,-p;p'^) f n±2^ \ 

9{p; p^)9{(jj; pYt] 4 ^ 

r/^(l-C) 

where we used Lemma [9. II and (14. 2p . Similarly, when n is odd, 

Xnip) = ^ ' 

^ * (1-C) 



Ynip) = -2 



7^^(1-0(1 + 0^ 

We now compare these expressions with (15.70 . For instance, when n = mod 6, 

Pn-l(0 - C^Pn-l(0 Tqn{T) 



-Xri 



4^(1-^2) r(p) 
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By Corollary 19.51 and (19 .4^ , this can be written as 

Pn-l(C) - C^Pn-l(C) = (1 - C)'(C' + 4C + 1)"'"^ 



4 

2 



n -| 

C(C + i)' \^~' /2(C2 + 4C + 1 



2(c2 + 4c + i)3); c(c + i 

Since, by Corollary 18 ■4[ p„_i(0) = 1, it follows that g„ is monic of degree n^/12 — 1, 
and that the final equation of Theorem 13.61 holds. Repeating the same argument 
for all cases, one obtains Theorem 13.21 and Theorem 13.61 

Finally, we comment on the deduction of Corollary 13.31 from Theorem 13.21 This 
is a tedious exercise, and we will only explain the case when n = 1 mod 6 and we 
consider squares of colour 2. In this case, 

^•in / \ / , n(n + 2) I tnt\ ifl^l ~t" ^n^2 ~l~ ^1 ^2 /™s 

Zl^'ihMM) = (W2)T K^g„(T) -2^i ^-^ —rn{T) 

V ^2 12 

= (W2) ^ — T 12 -2 T 12 )+■■ 

\ ^2 ^2 

the ellipsis denoting lower terms in T, hence also in Writing 

^3 



this can be simplified to 



^2 I 17 5n2+8n-13 „2 „ 

(totl)^^t2 (to+tl)~(to + i?) + --- • 

Thus, the leading power of ^2 is (Sn^ + 8n — 13)/12, and the coefficient of 

n2+4n + 7 , ,„ 5r 
f ^ f~ 



n2+4^i + 7 5n2+8n + ll , 5n2+8n-13 



IS 

'(^2 - 9)/4\ - 9)/4 

A; k-2 

Repeating the same argument for each colour and each residue class of n mod 6, 
one obtains Corollary 13.31 

9.3. Proof of Proposition [3741 In this Section, we prove that r2n+i has integer 
coefficients. We need the following fact, which is generalized in Corollary 19.81 

[ n(n + 2) ~\ 

Lemma 9.6. The polynomial 2 \- 2 \p^{^C/ 2) has integer coefficients. 
For the proof, we will use a determinant formula for p„. 
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Lemma 9.7. The polynomials pn are given by 



2"'n".iO -l)!=(l + 2C)[*l (1 + 



2j_[{n-l)£ 



X det 



F{x,y,2C + l] 
Gix,y) 



Proof. Let -> 2C + 1, yi C/(C + 2) and ^ = 2C + 1 in^D- The left-hand side 
can be expressed in terms of the polynomials pn using fl8.3bp . We now apply the 
well-known identity 

j-^ (^etl<iJ<nif{x^,yj)) _ 1 f d'+^-^f{a,b) 



Simphfying the prefactor using (18. 8p . we obtain the desired result. □ 

Proof of Lemma 19. 6[ In Lemma 19. 7[ replace a; by a; + 2^ -|- 1 and ?/by?/ + C/(C + 2) 
inside the determinant, and then replace C, by C/2. Note that 

F(a;-fC + l,y + C/(C + 4),C + l;C/2) ^ (C + 2)(C + 4)A(x, Q 

G(x + C + 1, 1/ + C/(C + 4); C/2) C) ' ^ • ^ 

where A, 5 G Z[a;,y,C] and 5(0, 0,C) = C^(C + 2)^. It follows that the Taylor 
expansion of (19. 6p at x = ?/ = has the form 

(C + 4) c,,{0 , , 



^'(^ + 2)0:;,(c^(c + 2)^)' 

where C^i G Using (18. 2p . Lemma 19771 can then be written 

[Ml _ (-l)("^')deti<,,,<. (a-i,,-i(C)) 

-'/^nl.S,/^; — ran r(„_i)2-| ■ 

^"^(^^ _^ 2)"("-i)(C + 1) (( + 4)"(""i)+L— 4— J 

The right-hand side is a quotient of two polynomials with integer coefficients. 
Since the denominator is monic, the left-hand side has integer coefficients. □ 

Proof of Proposition 13.41 Let n be even, and write (13. 4 p in the form 



n{n+l)/2 [n(n+2)/12] 



k=0 k=0 
[(n-2)(n+4)/12] 

- E |c'="^(C + ir^^(C^ + 4C + l)"^^-^^ (9.7) 

fc=0 

Here, 6^ are coefficients of the polynomial 2 2 pn{(/2), which are integer by 
Lemma [9^ The numbers Ck are coefficients of qn+i, which are also integers, see ^ 
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We need to show that dk are integers. To this end, let < m < [(n — 2)(n + 4)/12], 

ri(n + 2) 

pick out the coefficient of C ^ ™ on both sides of (19.71) . and multiply the result 
by 2™. This leads to a triangular system of the form 

fO, 0<m<f, 

Cm - dm + 2^{1^Cj + Zdj) = < n ^ ^ ^ r (n-2)(n+4) l 

j^m [On(^_^, 2 - - L 12 J- 

By induction on m, it follows that each dm is an integer. □ 

Proposition 13.41 implies an integrality result for the coefficients of p„. Although 
it seems to be far from sharp, it may still have some interest. 

Corollary 9.8. Let p^iC) = EfciV^^^^ "fcC"- Then, 2^'^ak G Z, where 

min ( mf-] , fc, ^ ~k), n odd. 



Proof. When n is even, this follows easily from (19. 7p . using that 2'^ ^ 2 ^ bk = 
and that Ck and dk are integers. For odd n, the proof is similar. □ 

10. Thermodynamic limit 

In this Section, we give a formal derivation of our conjectured expression (I3.12p 
for the free energy. Essentially, we assume that the free energy can be expressed in 
terms of a formal power series in the variable and then deduce (13 . 1 21) by formally 
taking — )■ 00 in Proposition 18.31 To turn this into a rigorous proof would require 
strong analyticity assumptions that seem difficult to verify a priori. However, one 
can observe that (I3.12p gives the correct result not only as C — ^ 0, but also at C = 1 
(because of (12.51) and (I3.14p ) and a.s C, ^ 00 (since it has the correct behaviour 
under ^ 1— )■ C~^)- It also seems to agree with numerical experiment. 

We will assume that 

log(Pn(C)) 1 /^N nni^ 

lim = log^(C) (10.1) 

n— ^00 n 

exists in the algebra of formal power series in C,. We also assume that the conver- 
gence is regular in the sense that if we factor 

Pn(C)=^7(CrVn(C), (10.2) 

then 

n^oo 0„+2(C)0n-2(C) 

still in the sense of formal power series. As a weak motivation for this assumption, 
we note that, by Proposition 18.91 and (I3.14p . (/)„(1) ~ Cn°'(3"', where C is different 
for even and odd n. Thus, (110. 3p holds pointwise a.t ( = 1, whereas the limit 



n— ^-oo 



0n+l(C)0n-l(C) 
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does not exist. 

Let us now apply Proposition 18.3^ which we write in the form 

(/-\ 1+X(" odd) 

where 0{("'^^) is a polynomial divisible by C"^^- Iterating this gives 

(1 + C)V+2(CK-2(C) = (1 + 20' (i + PniC? + 

Applying the factorization (110.21) we obtain 

(1 + 20^(1 + 1)' UC? , c?(C) 



^^^^ (1+C)« 0.+2(C)0n-2(C) ' (l + C)^^7(C)^">n+2(C)0n-2(C)^ 

Letting n — )■ oo and using (110. Sp we conclude that 



1 + C 

Replacing ( by in (110. ip then gives 

hm = log^(C) (10.4) 

?i--s>oo n 

Finally, using (110. ip and (110. 4p in Corollary 19.41 gives 

r log fa. ^1.^2) , (9{C)_\ 1, ,,x^i faiO 
hm = log = - log(totit2) + log 



n2 ^\7]y^e{X^;p3)J 3 ^^r7i/4 

By Lemma [9. 2 [ this may be written as in (I3.12p . 
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